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SECTIn^J  I 

TNTRODllCTinN 

This  report  discusses  testing  and  analysis  Arthur  D  Little,  Inc  fADL)  conducted 
under  the  Rotary  Reciprocating  Refrigerator  (RCURE)  Development  Program  Fll.  The 
RCUBE  cooler  operates  on  a  reverse-Brayton  cycle  to  provide  cooling  at  60  K  and  12 
K.  A  schematic  of  the  cycle  is  shown  in  Figure  1. 


Compressor  Assembly  — 

Heat  Exchanger  1 
Heat  Exchanger  W  2 

Heat  Exchanger  if3 
Expander  Assembly  _ 


Figure  1.  Schematic  of  two-expander  reverse-Brayton  cycle 

The  counterflow  heat  exchangers  are  vital  cooler  components  --  typically  each  is 
required  to  have  an  effectiveness  on  the  order  of  0.9R.  Because  of  the  large 
required  effectiveness,  the  heat  exchangers  make  up  a  significant  portion  of  the 
cooler  mass  and  volume.  In  order  to  reduce  heat  exchanger  mass  and  volume,  compact 
heat  exchangers  composed  of  perforated  plates  separated  hy  G-10  spacers  were 
developed.  Unfortunately,  the  organic  materials  in  the  heat  exchangers  can  outgas 
and  contaminate  the  working  fluid  which  reduces  system  lite. 
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Since  system  ii^e  is  a  primary  concern,  an  efi'ort  was  undertaken  to  develop 
all-metal  compact  heat  exchangers  to  replace  the  PCUBE  heat  exchangers  which 
contained  organic  materials.  The  heat  exchangers  would  he  composed  of  copper 
perforated  plates  separated  by  stainless  steel  e-seal  spacers.  They  would  still  be 
a  two-pass  counterflow  design  with  the  high  pressure  stream  in  the  circular  center 
section  and  the  low  pressure  stream  in  the  outer  annular  section  ^see  Figures  ’  and 
3). 

Note  that  the  e-seal  spacer  has  two  functions.  It  both  separates  the  two  ^low 
paths,  and  inhibits  longitudinal  conduction  between  ad.iacent  plates.  Because  of  the 
e-seal  geometry,  the  physical  plate  separation  is  0.087"  while  the  conduction  path 
lenqth  is  0.156"  fsee  Figure  4). 

High  Pressure  Low  Pressure 


Figure  ?.  Side  view  of  sample  heat-exchanger  with  six  plates  and  five 
e-seal  spacers 
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Low  Pressure  Stream 


Fiqure  3.  Top  view  of  plate  and  e-seal  spacer 


.087" 


Perforated  Plate 


'hole 


Figure  4.  E-seal  geometry  and  hole  misalignment 


Although  the  all-metal  heat  exchanger  reduces  contamination  and  therefore 
increases  system  life,  it  also  results  in  a  larger  heat  exchanger  for  a  given 
ef'^ectiveness  because  it  doesn't  provide  as  much  thermal  isolation  between  the 
plates.  Recause  system  life  is  considered  so  ’nportant,  the  increase  in  system  mass 
and  volume  is  acceptable. 
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RCIIBE  system  requirements  dictate  the  minimum  heat  exchanger  e'FFectiveness  and 
the  maximum  acceptable  pressure  drop  across  the  heat  exchangers.  Tn  order  to  design 
the  heat  exchangers,  heat-transfer  and  pressure  drop  data  had  to  be  gathered  for 
various  plate  geometries  and  flow  conditions.  The  required  effectiveness  o-f  the 
balanced,  counterflow  heat  exchangers  dictates  the  number  ot  heat-transfer  units 

for  the  entire  heat  exchanger.  Consequently,  N^^/plate  was  chosen  to  compare 
the  heat-transfer  characteristics  of  various  plate  geometries.  The  friction  factor 
was  chosen  to  compare  the  pressure  drop  characteristics.  Knowing  the  N^^/plate  and 
friction  factor,  the  number  of  plates  needed  to  achieve  the  required  effectiveness 
and  the  associated  heat  exchanger  pressure  drop  can  be  determined. 

Next,  test  methods  to  determine  N^^/p’ate  and  friction  factor  were  chosen.  An 
isothermal  pressure  drop  test  was  used  to  find  friction  factors  while  the 
single-blow  transient  test  method  [2]  was  used  to  gather  heat-transfer  data. 

The  single-blow  transient  test  method  allows  a  core  of  plates  separated  by 
e-seals  to  be  tested  using  only  one  flow  path  (see  Figure  F).  The  room  temperature 
core  is  placed  into  a  heated  air  flow  and  the  air  temperature  downstream  from  the 
core  is  recored  until  a  steady  state  is  reached.  The  core  can  be  determined 
from  the  maximum  slope  of  the  time-temperature  plot  and  the  longitudinal  heat 
conduction  in  the  core. 

Flow  Direction 


Figure  5.  Configuration  of  test  core 
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Since  each  plate  in  the  core  is  isothermal  at  any  given  time  during  this  test 
the  N^^/plate  data  reported  is  for  100^  fin  effectiveness.  In  order  to  use  the  data 
to  design  actual  counterflow  heat  exchangers,  the  "^in  effectiveness  must  be 
considered  f see  Appendix  C). 

Eleven  different  plate  geometries  were  tested  and  compared.  The  plate 
Geometries  that  provided  the  best  heat  transfer  and  friction  factors  were 
identified.  All  the  perforated  plate  test  sections  had  3"  diameter  copper  plates 
separated  by  stainless  steel  e-seal  spacers  and  were  designed  tor  a  helium  working 
fluid.  Two  test  sections  had  places  with  etched  holes  and  ?4.5^^  open  area),  one 

had  plates  with  punched  holes  f?3.6%  open  area),  and  eight  had  plates  manufactured 
using  electron  beam  drilling  (P0%  open  area). 

Originally  the  punched  and  etched  plates  were  tested  and  the  testing  indicated 
that  plates  with  smaller  hole  diameter  to  plate  thickness  (d/t)  ratio  had  better 
heat-transTer  characteristics.  Consequently,  ADL  searched  for  a  manufacturing 
method  which  could  produce  plates  with  smaller  d/t  ratios.  They  found  the  electron 
beam  drilling  process,  so  eight  more  cores  (0.154  1  d/t  i  0.750)  were  manufactured 
and  tested.  The  three  original  test  cores  had  40  plates  separated  by  39  e-seals. 

As  thicker  plates  were  manufactured,  the  cores  had  to  have  fewer  plates  to  fit  in 
the  switching  box  --  the  last  8  cores  had  24  to  32  plates.  Figure  6  shows  an 
assembled  test  section,  an  e-seal  spacer,  and  the  three  original  plates. 


Figure  6.  Assembled  test  section,  e-seal  spacer,  and  three  plate  types 
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Throiiahout  the  tests,  the  holes  in  the  test  cores  were  misaligned  (see  Figure 
4).  K  the  distance  between  the  plates  is  small  enough  the  alignment  nr 
misalignment  of  the  holes  is  expected  to  affect  both  heat-transfer  and  pressure  drop 
characteristics.  Misalioning  the  holes  is  expected  to  enhance  heat  transfer  while 
increasing  pressure  drop  and  aligning  the  holes  is  expected  to  minimize  pressure 
drop  at  the  expense  of  heat  transfer.  There  is  an  indication  fll  that  if  the  ratio 
of  plate  separation  to  plate  thickness  is  greater  than  one,  the  plates  won't  be 
influenced  by  adjacent  plates.  The  ratio  of  plate  separation  to  plate  thickness 
ranged  from  4.70  to  1.34  in  the  test  cores.  Consequently,  the  hole  misalignment 
shouldn't  affect  the  pressure  drop  through  the  cores.  Furthermore,  as  long  as  total 
exchanger  pressure  drop  remains  acceptably  small,  it's  desirable  to  enhance  heat 
transfer  at  the  expense  of  pressure  drop. 

This  report  describes  the  test  methods  and  apparatus  used  to  evaluate  and 
compare  the  plates.  Evaluations,  and  resulting  modifications,  of  the  test  method 
and  apparatus  are  also  discussed.  Data  gathered  using  the  final  method  were 
considered  accurate  enough  for  design  purposes  and  were  compared  to  theoretical 
results. 

Section  TT  discusses  the  characteristics  chosen  to  compare  the  test  sections 
and  specifies  the  measurements  necessary  to  calculate  the  plate  characteristics. 
Next,  this  section  describes  the  test  apparatus  and  procedure.  Finally,  the 
assumptions  of  the  single  blow  transient-test  and  the  test  conditions  that  closely 
approximate  these  assumptions  are  presented. 

Section  IIT  shows  the  poor  agreement  between  the  test  and  theoretical  results. 
The  reviow  of  tho  test  procedure  and  the  subsequent  modifications  o''^  the  test 
ipparatus  are  then  discussed.  Next,  Section  III  verifies  the  accuracy  of  results 
obtained  using  the  modified  test  apparatus.  Section  IV  then  presents  the  results. 
Sections  V  and  VI  cover  conclusions  and  recommendations. 

Appendix  A  presents  the  eguations  used  to  perform  data  reduction.  The  data, 
including  a  time-temperature  trace  from  the  modified  apparatus  and  results  are 
presented  here  and  Appendix  B  describes  the  accuracy  of  the  results. 
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The  theoretical  equations  used  to  check  the  experimental  data  are  derived  in 
Appendix  C.  Appendix  D  compares  the  steady  state  and  transient  analysis  to 
determine  the  insulation  requirements. 

The  theory  behind  the  single-blow  transient  test  method  is  clearly  presented  by 
Pucci  et  al  [2]. 

This  report  is  based  on  a  summary  report  [l]  documenting  tests  performed  by  Mr 
Richard  Hubbell.  The  tests  were  conducted  at  ADL  under  a  contract  to  develop  a 
two-stage  rotary  reciprocating  refrigerator.  Dr  Fowle,  a  consultant  to  ADL,  did  the 
necessary  analytical  work  to  support  the  tests. 
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SECTION  II 

TEST  SETUP 


1.  DATA 

Since  the  objective  oT  the  tests  was  to  compare  test  sections  composed  of 
different  plates,  a  method  of  comparison  had  to  be  chosen.  Since  the  number  of 
heat-transfer  units  per  plate  will  dictate  the  size  of  the  heat  exchanger  for 

any  required  heat  exchanger  effectiveness,  the  N^^/plate  was  selected  to  compare 
heat-transfer  data.  The  friction  factor  (f)  was  selected  for  the  pressure  drop 
comparison  and  data  were  gathered  for  Reynolds  numbers  ranging  from  4  to  666.  The 
Reynolds  number  used  for  all  calculations  ^hole  Reynolds  number)  was  calculated 
using  a  hole  diameter  as  the  representative  length,  and  the  mass  flow  rate  was  based 
on  the  minimum  free  flow  area.  The  following  properties  were  measured  and  used  to 
calculate  and  f:  the  mass  flow  rate,  the  pressure  drop  across  the  test  section, 
the  temperature  upstream  and  downstream  from  the  test  section  as  a  function  of  time. 

The  N^^  is  a  function  of  the  longitudinal  heat  conduction  (thermal  conduction 
in  the  solid,  parallel  to  the  flow  direction)  and  the  maximum  slope  of  the 
temperature  downstream  from  the  test  section  vs  time  plot.  The  relationship  between 
the  longitudinal  heat  conduction,  the  maximum  rate  of  change  of  temperature  with 
time  and  the  is  presented  in  tabular  and  graphic  form  by  Pucci  [2].  Note  that 
the  Colburn  modulus  can  be  calculated  given  the  heat-transfer  number  The 

Prandtl  number  (Np^),  and  the  plate  geometry  (A^/A)  using  the  following 
relationship: 


•i  =  (A^/Al 


fl) 


where 


j  =  Colburn  modulus. 

A  =  Minimum  free  flow  area, 
c 
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A  =  Total  heat-transfer  surface  arta. 
Nfu  =  Number  of  heat-transfer  units. 

Npj,  =  Prandtl  number. 


The  friction  factor  per  plate  can  be  calculated  given  the  mass  flow  rate  and 
pressure  drop  across  the  test  section  during  steady  state  operation. 


where 


AP  P 
_ c 

?  n  t  ih 


(2) 


f  =  Friction  factor. 

=  Hydraulic  diameter. 

AP  =  Pressure  change  across  the  test  section. 

p  =  Density  of  working  fluid. 

A  =  Minimum  free  flow  area, 
c 

n  =  Number  of  plates, 
t  =  Plate  thickness, 
tfi  =  Mass  flow  rate. 

Pc  =  Proportionality  factor  in  Newton's  second  law. 


h  c 

9  ^  is  a  physical  property  of  the  test  section,  and  the  temperature, 
p'r^essure  change,  and  volumetric  flow  rate  were  measured. 
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DESCRIPTION  OF  APPARATUS 


Switching  Box 


(Air  or  He) 


Figure  7.  Schematic  of  test  apparatus  (dummy  core  in  place) 


AFWAL-TR-88-3in2 


Figure  8.  Schematic  of  switch  box  operation 


The  test  apparatus  (Figure  consists  of  a  source  of  compressed  air  (or 
helium),  a  resistance  heater  followed  by  a  flow  straightening  section  and  the  test 
section.  The  test  section  is  in  a  switching  box  ^Figure  81  which  allows  it  to  be 
moved  quickly  into  the  stream  of  heated  fluid.  When  the  test  section  is  not  in 
the  system,  a  dummy  core,  which  causes  the  same  pressure  drop  as  the  test  section, 
is  in  place.  Temperature  and  pressure  readings  are  taken  downstream  from  the 
straightening  section,  and  the  volumetric  flow  rate  is  measured  upstream  from  the 
straightening  section  using  a  rotometer.  The  outlet  temperature  is  measured  after 
the  switching  box,  and  the  outlet  pressure  tap  is  far  enough  downstream  from  the 
test  section  ("^3  in)  to  allow  for  pressure  recovery. 

3.  TFST  PROCEnilRE 

The  test  procedure  used  to  collect  the  dafa  follows: 
a.  Number  of  Heat-Transfer  Units  ^N^^) 

(ll  Establish  a  uniform  temperature  in  the  test  core  by  blowing 
room-temperature  compressed  air,  or  helium,  through  the  test  apparatus  with  the 
test  section  in  place.  A  uniform  temperature  (±n.l°F)  is  necessary  through  the 
entire  test  core  volume  to  ensure  repeatable  results. 
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(2)  Move  the  dummy  section  into  the  flow  to  allow  the  fluid  to  bypass  the 
test  section. 

(3)  Use  the  resistance  heater  to  increase  the  fluid  temperature  about 

?0°F. 

(4)  Maintain  the  heated  airflow  until  all  the  apparatus  components  reach 
a  stable  temperature.  The  time- temperature  plot  should  vary  less  than  1.0°F/hr 
once  this  steady-state  condition  is  achieved. 

(5)  Next,  move  the  room-temperature  test  section  into  the  hot  air  stream, 
and  record  the  air  temperature  leaving  the  test  section  as  a  function  of  time. 

The  pressure  drop  across  the  test  section  and  the  dummy  section  must  be  within 
0,1"  HpO  of  each  other  to  maintain  a  constant  flow  rate  through  the  apparatus 
before  and  after  the  test  section  is  moved  into  the  stream. 

f6)  Repeat  the  procedure  for  various  mass  flow  rates. 

b.  Friction  Factor  ff) 

Record  the  pressure  drop  across  the  test  section  for  steady-state, 
room-temperature  fluid  as  a  function  of  the  flow  rate. 

4.  THEORETICAL  ASSUMPTIONS 

The  transient-test  technique  is  based  on  assumptions  which  the  test  method  and 
apparatus  must  insure.  The  assumptions  and  the  test  conditions  that  closely 
approximate  each  assumption  follow: 

a.  The  fluid  temperature  and  velocity  profiles  are  uniform  across  the 
cross-sectional  area  of  the  test  section.  The  flow  straighteners  insure  a  uniform 
velocity  profile,  and  the  time-temperature  plot  shows  when  a  steady-state,  uniform 
temperature  is  reached. 
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b.  The  fluid  flow  rate  is  steady.  This  is  assured  by  keeping  a  large 
reservoir  of  fluid  and  insuring  that  the  pressure  drop  is  the  same  across  the  test 
and  dummy  sections. 

c.  Longitudinal  heat  conduction  is  minimal.  This  is  a  design  characteristic 
of  perforated-plate  heat  exchangers. 

d.  Heat  conduction  across  the  core  holder  is  negligible.  An  insulated  balsa 
wood  holder,  which  has  low  thermal  canacity  and  conductivity,  meets  this 
requirement. 

e.  The  matrix  (core)  heat-transfer  coefficient  is  independent  of  position 
and  time.  The  plate  construction  is  uniform,  and  heat-trans'^er  and  fluid  flow 
effects  are  distributed  uniformly  through  the  core. 

f.  The  thermal  properties  of  the  test  fluid  are  independent  of  temperature. 

A  20°F  temperature  change  assures  that  the  properties  do  not  change  substantially. 
This  condition  was  tested  by  comparing  the  N^^/plate  for  20,  30,  and  40°F  changes 
in  air.  The  results  were  within  4^  of  each  other  for  the  same  Reynolds  number. 

q.  The  gas  acts  as  a  perfect  gas.  At  ambient  temperature  and  pressure  (test 
conditions)  air  and  helium  can  be  considered  perfect  gases.  Also,  using  air  as  a 
test  gas  to  gather  data  on  helium  heat  exchangers  is  appropriate  because  the 
Prandtl  numbers  of  air  and  helium  are  very  close  (Np^(air)  =  0.72,  Np^(He)  = 

h.  The  thermal  capacity  of  the  gas  in  the  matrix  is  small  compared  to  the 
thermal  capacity  o^  the  matrix.  A  core  of  24  -  40  plates  was  used  for  the  test 
assuring  good  thermal  capacity. 

i.  Heat  conduction  in  the  gas  is  negligible  in  the  flow  direction.  The 
convective  heat  transfer  is  much  greater  than  any  gas  conduction  effects. 
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SECTION  ITT 

EVALUATION  OF  TEST  APPARATUS 


1.  COMPARISON  WITH  THEORY 


The  N  /plate  can  be  detennined  as  a  function  of  plate  geometry  and  hole 


tu 

Reynolds  number. 


The  equations  used  to  calculate  the  N^^/plate  assume  that 


1  <  Np^  ^  1000 


N 


Pr 


0.7 


and  the  section  has  a  uniform  wall  temperature.  Two  laminar  Flow  cases  were 
considered.  One  case  assumed  a  developing  velocity  profile  in  each  hole  and  the 
other  assumed  a  fully  developed  parabolic  velocity  profile.  The  results  of  the 
calculations  in  Appendix  C  follow. 

FOR  A  DEVELOPING  VELOCITY  PROFILE 


N,^/PLATE  = 


-0.232  0.520  0.121 

{l-a)a  (2.554+2.011  N 


ll 


Ke 


)  20.9 

—  +  - ST 


0.417 


Re 


(3) 


Re 


l+0.012[(d/t)Npgl 


0.8 


FOR  A  PARABOLIC  VELOCITY  PROFILE 


N^^/PLATE 


fl-a)a'^-^^2(2. 554+2. 011  0.269 


'Re 


'Re 


l+0.032r(d/tlNn 

(4) 


where 

=  Heat-transfer  number. 

0  =  Ratio  of  plate  open  area  to  frontal  area, 
t  =  Plate  thickness, 
d  =  Hole  diameter. 

Nro  =  Hole  Reynolds  number. 


13 


AFUAL-TR.R8-3102 


Comparing  the  predicted  values  with  the  test  results  (Figure  9)  showed  that 
the  predicted  N^^^/plate  are  ?.3  to  A. 7  times  the  test  results.  The  actual  and 
predicted  times  to  maximum  slope  ot  the  temperature  vs  time  curve  were  also 
compared.  The  time  to  maximum  slope  (0*)  is  based  on  the  mass  tinw  rate  and  the 
Fluid  and  heat  exchanger  properties. 


Where : 


m  Cp(air) 

It  =  Mass  Flow  rate  of  the  Fluid. 

n  =  Mass  oF  plates, 
m 

Cjjj  =  Specific  heat  of  plates. 

=  Specific  heat  of  air  (or  helium), 
pair  ' 


(B) 


Table  1  compares  the  calculated  and  measured  time  to  maximum  slope  for  runs  3 
and  S. 


0  (sec) 

RUN  ACTUAL  PREDICTED 
3  169  151 

5  87.5  74.1 

Table  1.  Comparison  of  calculated  and  measured  time  to  the  maximum 
cooling  rate 


The  actual  values  were  read  from  the  time  temperature  traces.  The  difference 
in  the  values  indicates  a  fluid  loss;  however,  it  is  very  difficult  to  determine 
the  time  to  maximum  slope  accurately  from  the  traces. 

Since  there  were  large  differences  between  the  actual  and  predicted  values  of 
both  the  N^^^/plate  and  time  to  maximum  slope,  the  test  method  and  apparatus  was 
suspect. 
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?.  REVIEW  OF  TEST  PROCEDURE 

Because  of  the  discrepancy  between  predicted  and  measured  values  of  N^^/plate 
vs  ''Re  and  time  to  maximum  slope,  the  test  procedure  was  reviewed.  The  followina 
aspects  of  the  test  were  checked: 

ASSUMPTIONS.  V/ere  the  assumptions  used  to  model  the  heat  exchanger  and 
predict  the  results  met  during  the  test’ 

INSTRUMENTATION.  Was  the  instrumentation  accurate  enough  to  provide  reliable 
data’ 

CALCULATIONS.  Were  the  methods  and  calculations  used  to  manipulate  the  data 
correct? 

The  instrumentation  and  calculations  were  fine;  however,  there  was  a  large 
apparent  heat  leak  that  could  have  been  due  to  either  heat  leaks,  fluid  leaks,  or 
some  combination  of  the  two.  The  test  apparatus  model  assumed  that  mass  flow 
remained  constant  (no  flow  leakage)  and  that  heat  conduction  across  the  core 
holder  is  negligible  (no  heat  leakage).  The  assumptions  upon  which  the 
transient-test  technique  is  based  were  not  met. 

a.  Heat  Leaks 

A  heat  leak  was  suspected  when  a  temperature  difference  was  noticed  across 

the  core  once  the  system  reached  steady  state.  The  instrumentation  was  checked  -- 

it  provided  accurate  sensing.  The  apparent  heat  leak  was  then  estimated  based  on 
the  following  core  information: 

T^.  =  P8°F  =  Inlet  temperature. 

T^  =  95.3°F  =  Outlet  temperature. 

’'AT  =  20°F  =  Initial  temperature  step  applied  to  the 

core. 

It  =  0.2?  Ibm/min  =  Mass  flow  rate. 

'A9  =  340  s  =  Time  to  reach  steady  state. 
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m  =  l.?87  Ibn  =  Mass  of  plates, 
s 

=  0.09  RTD/lbm  R  =  Specific  heat  of  plate. 

Cp(air)  =  0.?A  BTU/lbm  R  =  Specific  heat  of  air. 

If  the  temperature  difference  was  due  only  to  a  heat  leak,  the  heat  loss  rate 
can  be  found  from  the  following  relationship: 

0^  =  111  Cpfair)  (T.  -  T^) 

=  (0.??  Ibm/min)  (0.?4  BTH/lhm  R)  (98  -  95.3)F  ^eOmin/hr) 

=  8.B5  RTU/hr 

To  find  the  percent  heat  leak,  the  rate  of  heat  storage  must  be  estimated. 

The  heat  absorbed  by  the  heat  exchanger  from  the  start  of  the  test  to  steady  state 
(Q)  is  calculated  using  the  following  equation: 

0  =  m^c..  AT 
s  s 

=  (1.287  Ibm)  (0.09  BTU/lbm  R)  (20°F) 

=  2.316  BTU 

The  average  rate  of  heat  storage  (q^)  can  now  be  approximated: 

.  0  2.316  BTU 

9s  =  _  =  -  (3600  s/hr) 

A0  340  s 
=  24.52  BTU/hr 

The  resulting  apparent  heat  loss  rate  (Q^  =  8.55  BTU/hr)  is  35^1!  of  the  heat 
storage  rate  =  24.55  RTli/hr)  —  this  is  an  unacceptable  loss. 
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h.  Fluid  Leaks 

The  apparent  heat  leak  may  be  partially  due  to  fluid  leaks.  A  reasonable 
estimate  o-f  the  heat  leak  is  4.5  RTlI/hr.  If  the  heat  leak  was  actually  only  4.5 
BTU/hr,  the  fluid  leak  can  he  estimated  by  performing  an  energy  balance  on  the 
system. 


Figure  10.  Test  section  control  volume 

The  ^luid  loss  at  Section  1  is  the  critical  loss;  any  loss  at  Section  ?. 
should  not  significantly  a'Ffect  the  test. 

dE 

/g^dt  =/—  dt  +  m2^)CpT,dt  +  +  nij)CpTjdt 


where 

ifij  =  Mass  flow  into  test  section, 
ihj^l  =  Mass  flow  leaking  out  of  test  section  at 
the  entrance. 

iTtp  =  Mass  flow  out  of  the  test  section. 

liipi  =  Mass  flow  leaking  out  of  the  test  section  at 
the  exit. 

t  =  Time. 

E  =  Internal  energy. 

AiTii  =  trij^  finn%l  =  Percent  fluid  loss. 
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The  fnllowing  relatimship  results  from  the  continuity  equation: 


Now,  assuming  all  the  flows  are  constant,  for  the  steady-state  condition  the 
percent  fluid  loss  can  he  found  directly  from  the  following  relationship: 

^1 

=  1 - (lOO^K) 

ni,Cp(Ti-T^) 

j  _  _ 4.5  BTIJ/hr _  hr 

(0.??  Ibm/min)  (0.24  BTU/lbm  R)(2.7F)  60  min 

=  47% 

The  calculations  indicated  the  need  to  modify  the  apparatus  to  reduce  both  heat 
and  fluid  leaks. 

3.  MOOTFTCATTON  OF  TEST  APPARATUS 


The  test  apparatus  was  modified  using  a  seal  tight  design  to  reduce  fluid  leaks 
and  more  effective  insulation  to  reduce  heat  leaks.  The  transient  behavior  ot  the 
insulation,  not  the  steady  state  behavior,  determined  how  much  insulation  was 
needed.  Figure  8  compares  the  estimated  total  heat  loss  vs  t,.Tie  calculated  using 
steady  state  and  transient  analysis.  Measured  values  of  the  apparent  heat  leak  are 
also  presented.  Runs  25  and  30  were  performed  using  the  modified  apparatus.  Refer** 
the  modification,  the  discrepancy  between  the  measured  and  predicted  values  (using 
both  steady  state  and  transient  analysis)  of  the  heat  leak  were  large.  After  the 
modification  the  heat  leak  for  long  runs  had  decreased  considerably  and  the  measure  ' 
and  predicted  values  fusing  transient  analysis)  were  much  closer. 
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Elapsed  Time 't  (sec) 

Figure  11.  Comparison  of  measured  and  predicted  heat  leaks 

The  modified  apparatus  reduced  heat  and  flow  leaks,  producing  results  accurate 
enough  for  design  evaluations.  The  accuracy  of  the  apparatus  can  be  demonstrated  by 
comparing  the  actual  and  theoretical  apparent  heat  leaks  and  time  to  the  maximum 
temperature  slope  for  various  runs. 

The  heat  loss  from  the  system  was  found  using  an  energy  balance  assuming  no 
fluid  leakage. 


Figure  1?.  Test  section  control  volume 
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where 


til 

n,  =  /Qdt  =  dt  +  dt  -  /  ihc^T.dt 

loss  n  ndt  ft  P  ?  ft  P  1 


=  ■’’’mF-m  ftx  -  t.)  +  ltlC_  /  (T,  -  T  )dt 
mtn  T  1  1  c 


=  Mass  of  core, 
m 

c„  =  Specific  heat  of  core, 
m 


t^  =  Temperature  of  core  at  final  steady  state, 
t^.  =  Temperature  of  core  before  heatinq. 

-  Tj,)dt  =  Area  under  the  temperature  vs  time  plot, 

T  =  Time  for  the  run  to  reach  steady  state. 


The  test  generates  data  to  calculate  the  actual  heat  leak.  The  predicted  heat 
loss  for  infinite  insulation  can  be  estimated  using  tabulated  values  of  heat  loss  as 
a  function  of  insulation  properties  and  the  time  needed  to  reach  the  steady-state 
condition.  The  actual  time  to  maximum  slope  can  be  read  from  the  time-temperature 
plot  and  equation  (5)  generates  the  predicted  time  to  maximum  slope. 


The  results  of  the  comparisons  for  test  runs  ?5  and  30  are  presented  in  Table  ?. 
The  runs  show  excellent  agreement  between  the  predicted  and  measured  values  and 
indicate  minimal  heat  and  fluid  losses. 


RUN 

T  (sec) 

0  (RTII) 

Measured  Predicted 

6  (sec) 

Measured  Predicted 

?5 

117 

0,25 

0.17 

38.0 

38.9 

30 

420 

0.45 

0.35 

293 

298 

Table  2.  Comparison  of  actual  and  predicted  apparent  heat 
leak  and  time  to  maximum  slope 


21 


AFWAL-TR-R8-310? 


SECTION  IV 
TEST  RESULTS 

The  three  heat  exchanger  core  types  f punched,  etched,  and  electron  beam  drilled 
^EBnl)  were  tested  Tour  ways.  The  punched  and  EBD  plates  have  a  breakout  caused  by 
manufacturing  (See  Figure  13).  Tests  were  run  with  the  breakout  facing  both 
upstream  and  downstream.  The  orientation  of  the  core  affected  the  N^^/plate  but  did 
not  affect  the  friction  factor.  Since  the  etched  plates  do  not  have  breakouts, 
orientation  was  not  considered.  The  etched  plates  were  tested  with  both  helium  and 
air  to  insure  there  was  not  a  significant  Prandtl  number  effect — there  wasn't. 

Figure  14  presents  both  the  theoretical  and  measured  friction  factors  as  a 
function  of  hole  Reynolds  number.  The  theoretical  results  were  calculated  using 
equation  (C44)  and  are  presented  in  Tables  Cl  and  CP.  The  theoretical  results  show 
that  the  friction  factor  decreases  as  the  percent  open  area  increases  and  as  the 
hole  diameter  to  plate  thickness  ratio  decreases.  The  experimental  results  indicate 
that  the  friction  factor  is  more  sensitive  to  the  hole  diameter  to  plate  thickness 
ratio  than  the  percent  open  area.  As  the  hole  diameter  to  plate  thickness  ratio 
decreases,  the  friction  factor  also  decreases. 

Fiqure  15  shows  that  the  friction  factor  data  didn't  depend  on  the  test  gas 
(air  nr  helium)  or  the  core  orientation. 

Figure  16  presents  the  theoretical  and  measured  heat-transfer  data.  The 
theoretical  results  presented  assume  a  parabolic  veloci^v  profile  in  the  holes,  were 
calculated  using  equations  in  appendix  C,  and  are  tabulated  in  Tables  C3  and  C4. 

The  theoretical  results  indicate  that  N^^/plate  increases  with  decreasing  d/t  ratio. 
The  experimental  results  indicate  that  for  hole  Reynolds  numbers  above  30,  the 
Nty/plate  increases  as  d/t  decreases.  Also,  the  slope  of  the  N^^/plate  vs 
decreases  as  d/t  increases. 

Figure  17  demonstrates  the  dependence  of  the  heat-transfer  data  on  plate 
orientation  (or  hole  shape).  Both  the  punched  and  EBD  plates  were  oriented  in  the 
core  with  all  the  hole  breakouts  facing  the  same  direction.  When  the  core  was 
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Figure  13a.  Photograph  of  punched  and  etched  plate  cross  section 
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Figure  13b.  Photo^aph  of  EBD  plate  cross  section  (t  =  .065",  d  =  .010") 


Comparison  of  theoretical  and  actual  friction  factor  data 
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Figure  15.  Friction  factor  as  a  function  of  Reynolds  number 
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theoretical  and  actual  heat -transfer  character! 
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Figure  17.  Effect  of  plate  orientation  on  Ntu/plate 
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oriented  so  the  fluid  in  the  holes  followed  a  converging  path,  the  N^^/plate  for  a 
given  was  higher  than  for  a  core  oriented  so  the  fluid  followed  a  diverging 
path.  This  small,  but  measurable,  increase  in  N^^/plate  is  probably  due  to  less 
flow  separation  in  the  holes. 

The  heat-transfer  data  didn't  depend  on  the  test  gas  (Figure  18). 

Note  that  the  theoretical  heat-transfer  predictions  are  not  accurate.  Both  the 
predicted  slope  and  intercepts  are  off  resulting  in  predictions  that  are  too  high 
(by  as  much  as  a  factor  of  ?0  for  low  and  low  d/t).  Also,  the  model  failed  to 
predict  that  larger  d/t  ratios  result  in  larger  drops  in  N^^/plate  as  Np^  increases. 
The  discrepancies  are  probably  due  to  uncertainty  over  actual  fluid  flow  in  the 
holes.  Clearly  hole  shape  affects  the  heat-transfer  data  and  the  model  assumes 
cylindrical  holes  with  no  flow  separation.  Actually  the  holes  are  not 
cylindrical  and  there  may  be  flow  separation  which  accounts  for  the 
lower  measured  N^^/plate  values. 


Figure  18.  Comparison  of  helium  and  air  heat- transfer  data 
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SECTION  V 

CONCLUSIONS 

1.  The  modified  apparatus  using  the  single  blow  transient  test  method  to  determine 
the  heat-transfer  characteristics  of  perforated-plate  heat  exchangers  produces 
results  accurate  enough  for  design  evaluations. 

2.  The  transient  test  method  is  somewhat  insensitive  to  heat  leaks  when  used  to 
determine  the  average  heat-transfer  characteristics  of  perforated-plate  heat 
exchangers.  The  test  technique  is;  however,  very  sensitive  to  heat  leaks  (and  all 
assumptions)  on  an  individual  basis.  If  the  assumptions  aren't  insured  by  the  test 
method,  the  individual  results  aren't  repeatable. 

3.  Transient  behavior  dictates  the  amount  of  insulation  needed  for  the  test 
section,  not  steady-state  conditions. 

4.  Air  can  be  used  to  test  heat  exchangers  designed  for  helium  systems. 

5.  Hole  shape  affects  the  heat-transfer  characteristics.  The  increase  in  N^^/plate 
when  the  punched  and  ERD  plates  were  tested  with  the  breakout  facing  upstream 
indicates  that  a  converging-diverging  hole  may  provide  better  heat-transfer 
characteristics. 

6.  Heat-transfer  characteristics  depend  strongly  on  the  hole  diameter  to  plate 
thickness  ratio.  For  hole  Reynolds  numbers  greater  than  30,  the  N^^/plate  increases 
as  d/t  decreases.  Also,  the  N^^/plate  doesn't  decrease  as  quickly  with  increasing 
Np^  for  smaller  d/t  ratios. 

7.  The  model  used  to  determine  the  N^^/plate  doesn't  completely  describe  the 
system.  The  correlation  between  the  theoretical  and  experimental  results  is  poor, 
but  is  better  assuming  a  parabolic  velocity  profile  instead  of  a  developing  velocity 
profile  in  the  holes.  The  difference  between  the  predicted  and  actual  values  is 
probably  due  to  uncertainty  over  actual  flow  conditions  in  the  holes. 
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SECTION  VT 

RECOMMENDATIONS 

1.  Further  investigation  should  be  conducted  to  determine  the  effect  of  hole  shape 
on  N^j^/plate.  Also,  manufacturing  techniques  should  be  investigated  to  find  ways  to 
manufacture  plates  with  a  smaller  d/t  ratio. 

?.  Use  an  analog  to  digital  converter  and  a  minicomputer  to  record  the  temperature 
downstream  from  the  test  section.  This  will  reduce  the  possible  error  in  reading 
the  slope  and  will  reduce  the  corresponding  error  in  the  N^^/plate  vs  hole  Reynolds 
number  from  ±  6%  How  Np^)  or  ±  3%  (higher  Np^)  to  less  than  ±  \%. 
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APPENDIX  A  -  DATA  REDUCTION 


The  test  sections  that  were  tested  are  shown  in  Figure  1.  Temperature, 
pressure,  volumetric  flowrate,  and  the  time-temperature  history  were  recorded  for 
each  run.  This  data  was  reduced  using  the  following  method  to  calculate  the  Colburn 
modulus  and  friction  factor  vs  hole  Reynolds  number.  The  calculations  are  presented 
in  the  ■'"ol lowing  order:  mass  flow,  hole  Reynolds  number,  number  nf  heat-transfer 
units,  Colburn  modulus,  and  friction  factor. 


1 .  MASS  FLOW 

Start  with  the  definition  of  mass  flow: 


(Al) 


(^?) 
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Where : 

V  =  Velocity  of  fluid  entering  test  section. 
=  Open  area  of  plate. 

=  Mass  flow  rate. 

Vx  =  VA  =  Volumetric  ^low  rate, 
f  c 

R  =  Gas  constant. 

I  p  =  Density. 

P  =  Pressure. 

T  =  Temperature. 


The  volumetric  flow  rate  is  available  as  a  percentage  of  a  full-scale 
reading. 


Let: 

Gp  =  Flow  meter  full-scale  flow. 

Sp  =  Scale  reading. 

R«,„  =  53.34  ft  Ibf 

TBiT" 

Rm  =  386.0  ft  Ibf 
TbiO" 


Vp 
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For  air: 

(*i  (Ibm/min)  =  1.6444 


/  ^Ibm^  ^Ibf ^  I 

Gp(ft3/niin)  fA4) 

T  (R) 

NOTE:  Use  upstream  fluid  properties. 

2.  HOLE  REYNOLDS  NUMBER 


For  hel ium: 

it  flbrn/min)  =  0.611  S^ 


(A3) 


Again,  start  with  the  definition  of  the  hole  Reynolds  number: 


’Re 


u 

m  D. 


A^y 


(A5) 
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FROM  THE  OEFINTTION  OF  A  : 

c 

_  a(^DO 

■  4 

,  ,  9 

=  aiT(  ^  )  in^ 


P’oTry 


4rn  nhm/rm'n)  0^  (in) 

OTTfOin^)  y  Ohm) 
sec  ft 


min  l?in 

60sec  ft 


?.829  X  10“^  iti  (Ibm/min)  (in) 

.0  y  (1bm/s  ft) 


Where: 

ih  =  Mass  flow  rate. 

=  Hydraulic  diameter  of  hole. 

n  =  Plate  diameter. 

d  =  Hole  diameter. 

a  =  Ratio  of  the  plate  open  area  to  frontal  area, 
y  =  Fluid  viscosity. 

The  Reynolds  number  is  based  on  the  charateristics  of  one  hole.  Note  that 
the  hydraulic  diameter  of  any  cylindrical  pipe  is  the  pipe's  diameter. 
Consequently,  the  equation  for  Reynolds  number  can  be  rewritten  as  follows. 


_2 

?.8?q  X  10  m( Ibm/min)  d(in) 
a  y  Mbm/s  ft) 


(A61 
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3.  NUMRFR  OF  HFAT-TRANSFER  UNITS 

Pucci  [2]  presents  the  N^^^  as  a  function  of  the  maximum  slope  and 
longitudinal  conduction  parameters  both  graphically  (Figure  3[2])  and  tabularly 
(Table  1[2]).  To  find  the  the  maximum  slope  and  longitudinal  conduction 
parameters  must  be  calculated. 

a.  MAXIMUM  SLOPE  PARAMETER 

The  maximum  slope  parameter  is  defined  as  follows; 

dt^j,/de  I  max 

MAX.  SLOPE  PAR  =  - ^ -  (A7l 

ffi  ^pAt^^ 

Where : 

The  maximum  slope  of  the  trace  of  outlet 
temperature  vs  time. 

Step  change  imposed  on  inlet  gas  temperature. 

M^  =  The  mass  of  the  heat-transfer  core. 

=  The  specific  heat  of  the  heat-transfer  core. 

Cp  =  The  specific  heat  of  the  fluid. 

It)  =  Mass  flow  rate  of  the  fluid. 

The  mass  of  the  core  can  be  determined  from  the  plate  geometry  and  the  number 
of  plates. 

M  =  (l2^Ul-a)tnp  (A8) 

s  4 
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where 

n  =  Plate  diameter. 

a  =  Ratio  of  plate  open  area  to  frontal  area, 
p  =  Density  of  plate, 
t  =  Plate  thickness, 
n  =  number  of  plates  in  the  core. 

b.  Lonqitudinal  Conduction  Parameter 


k|^  =  9  RTIJ/hr  ft  R  fstainless  steel) 

=  0.015  in2 
L|^  =  4.3  in 


Figure  Al.  E-seal  dimensions 
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Wal  =  '  BTU/hrftR 


(Inconel ) 


esea 


^  ro.n08’')  tt 


'"eseal  ""  (0.156"  in/e-seal)  (39  e-seal) 


=  ^9  BTII/hr  ft  R)  (n.OlSinM  ft  hr 

4.2in  12in  60min 


=  4.464x10“'’  BTU/tnin  R 


*^eseal  %seal 


"eseal 


^  (7  BTU/hr  ft  R)  (2.25in)  (O.OOR)Tr  ft  hr 

(0.156  1n/e-seal)  (39  e-seal)  12in  60min 

=  9.036  X  10"^  BTU/rnin  R 


(4.464  +  9.036)  x  10“^  BTU/min  R 
It  (Ibm/min)  (0.24  BTU/lbm  R) 

5.63  X  10“^  (For  air) 


(For  helium) 


It  ( 1  bm/mi  n ) 

1.09  X  10“^ 
It  (Ibm/min) 


4.  COLBURN  MODULUS  (.1-factor) 


The  definition  of  the  Colburn  modulus  follows: 


h _ 

0  c_ 


(AlO) 
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Rut  N 

tu 


M. 

Itic„ 

p 


^  _  ^tu  111 
A 


.i 


m  2 '3 

N 

G  c  A  ”Pr 
P 


2/3 


rAii) 


Where: 

^u 

=  Heat-transfer  units. 

Ac 

=  Minimum  free  flow  area. 

A 

=  Total  heat-transfer  surface  area. 

=  Prandtl  number. 

j 

=  Colburn  modulus. 

5.  FRICTION  FACTOR 


nt?pV2 

But  iTi  =  PA^V 

2nt  it|2 


(A12) 
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Where: 

f  =  Friction  factor. 

,  AP  =  Pressure  change  across  the  test  section, 
n  =  Number  of  plates, 
t  =  Plate  thickness. 

~  Proportional ity  factor  in  Newton's  second  law. 

Tables  A1-A5  present  the  reduced  data.  Tables  A1  and  A2  present  the  friction 
factor  data.  Tables  A3  and  A4  present  the  heat-transfer  data  collected  using  the 
modified  apparatus,  and  Table  A5  presents  the  heat-transfer  data  collected  using 
the  unmodified  apparatus.  The  EBD  and  punched  core  tests  were  run  with  different 
plate  orientations.  The  core  descriptions  in  Tables  A1-A5  refer  to  the  following 
plate  orientations : 


A1 ternating 


H  i 

iS 


Enhanced 

Breakout  facing  upstream 


ss 

i  i 


ii 
$  i 
i  i 

Figure  A1  is  a  sample  time/temperature  plot  for  one  run. 


Not  Enhanced 

Breakout  facing  downstream 
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■ 

/Ibmx 
iti  min  ’ 

. 

Apf'H^O) 

NRe 

f 

GAS 

Core 

Description 

1 

0.09 

4.29 

9.52 

Air 

Punched 

2 

0.049 

0.13 

8.57 

3.44 

Plate 

3 

0.19 

76.0 

17.15 

1.26 

4 

0.173 

0.38 

76.0 

29.99 

0.82 

Breakout 

5 

0.247 

0.59 

76.0 

42.85 

0.63 

■facing 

upstream 

6 

0.119 

0.24 

76.0 

20.57 

1.10 

d=0.081" 

0.237 

0.57 

76.0 

41.13 

0.66 

t=0.021" 

0.356 

1.03 

76.0 

61.70 

0.53 

^0=0.236 

0.474 

1.70 

76.0 

82.25 

0.49 

* 

D 

0.911 

4.65 

76.0 

158.11 

0.36 

n 

0.075 

H 

75.0 

13.0 

1.39 

0.100 

75.0 

17.34 

1.04 

0.120 

75.0 

1.04 

Breakout 

14 

0.239 

75.0 

41.61 

0.62 

■facing 

15 

0.335 

0.91 

75.0 

58.24 

0.52 

downstream 

16 

0.479 

1.65 

75.0 

83.20 

0.47 

BMl 

mEsmmm 

19 

0.018 

0.77 

70.0 

1.43 

Plate 

20 

0.072 

1.75 

70.0 

0.81 

21 

BBSS 

3.20 

70.0 

76.80 

0.66 

d=0.0285'' 

2? 

4.30 

70.0 

89.60 

0.65 

t=0.0185" 

23 

■H 

6.50 

70.0 

115.21 

0.60 

=0.09 

24 

0.025 

18.04 

i 

25 

0.050 

36.06 

26 

0.075 

75.0 

54.12 

0.77 

27 

0.100 

75.0 

72.16 

0.68 

28 

0.125 

75.0 

90.27 

0.61 

29 

0.150 

0.73 

75.0 

108.31 

0.56 

30 

0.175 

0.93 

75.0 

126.35 

0.52 

31 

0.183 

0.99 

75.0 

131.76 

0.51 

32 

0.192 

1.16 

75.0 

138.62 

0.55 

33 

0.240 

1.70 

75.0 

173.25 

0.51 

34 

0.360 

75.0 

259.92 

0.48 

35 

0.432 

75.0 

311.87 

0.47 

36 

0.480 

6.31 

75.0 

346.51 

0.48 

37 

0.923 

18.20 

75.0 

666.1 

0.37 

_ 

_ 

Table  Al.  Friction  factor  data  -  punched  and  etched  plates  (40  plates/cnrel 
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u 

— 

Vluidf°F) 

"r. 

f 

-  1 

GAS 

Core 

■ 

imn 

HI 

Description 

■MH 

msgm 

Etched 

1 

0.0508 

16.1 

2.727 

Plate 

2 

0.1016 

mSSm 

BSH 

32.21 

1.364 

3 

0.1523 

mSSM 

48.31 

1.03Q 

d=0.0295" 

4 

0.2031 

BS9 

69.0 

64.41 

0.828 

t=0.0185" 

5 

0.2539 

69.0 

80.52 

0.701 

a =0.245 

6 

0.0487 

KK9 

69.0 

15.46 

2.961 

7 

0.0975 

0.14 

69.0 

30.91 

1.480 

8 

0.1949 

0.33 

69.0 

61.82 

0.872 

q 

0.2924 

0.56 

69.0 

92.73 

0.658 

10 

0.3898 

0.85 

69.0 

123.60 

0.562 

11 

0.4873 

1.21 

69.0 

154.60 

0.512 

■ 

Table  Al.  Concluded 
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Run 

Hj 

APC'H^O) 

Vluid(°F) 

"Re 

f 

1 

0.929 

6.14 

mvstgm 

■Mia 

? 

0.742 

4.12 

70.0 

3 

0.481 

2.67 

70.0 

Biia 

4 

2.3 

70.0 

54.6 

BwB 

5 

0.385 

1.94 

70.0 

48.5 

6 

1.63 

70.0 

42.5 

7 

0.288 

1.34 

70.0 

36.4 

Boil 

8 

0.240 

1.07 

70.0 

30.3 

BBi 

q 

0.192 

0.82 

70.0 

24.2 

WeeM 

in 

0.144 

0.60 

70.0 

18.2 

n.3’8 

11 

0.096 

0.39 

70.0 

12.1 

0.480 

12 

0.21 

70.0 

6.1 

1.03 

13 

0.50 

70.0 

15.8 

0.363 

14 

0.100 

0.04 

70.0 

12.6 

0.442 

15 

0.075 

0.30 

70.0 

9.5 

0.605 

16 

0.050 

0.20 

70.0 

6.3 

0.908 

17 

0.025 

0.12 

70.0 

3.2 

2.2 

69 

0.^40 

llS.l 

6.121 

70 

0.749 

BIH 

94.2 

0.126 

71 

0.486 

61.1 

0.199 

72 

masm 

73.0 

55.0 

0.213 

73 

iiil 

2.74 

73.0 

48.8 

0.229 

74 

0.340 

2.30 

73.0 

42.7 

0.251 

75 

0.291 

1.88 

73.0 

36.6 

0.280 

76 

0.242 

1.50 

73.0 

30.5 

0.322 

77 

0.194 

1.16 

73.0 

24.4 

0.389 

78 

0.145 

0.83 

73.0 

18.3 

0.495 

79 

0.097 

0.55 

73.0 

12.2 

0.738 

80 

0.29 

73.0 

6.1 

1.557 

81 

0.70 

73.0 

15.9 

0.554 

82 

0.101 

0.55 

73.0 

12.7 

0.680 

83 

0.40 

73.0 

9.5 

0.879 

84 

0.28 

73.0 

6.3 

1.384 

85 

mi 

0.14 

73.0 

3.2 

2.770 

Table  Friction  factor  data  -  ERD  plates 
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Run 

n 

APfH^O) 

"Re 

f 

GAS 

Core 

Description 

5? 

0.837 

8.02 

157.3 

0.176 

Air 

27  Plates 

53 

0.742 

5.89 

139.4 

0.164 

54 

3.68 

90.4 

0.243 

d=0.015” 

55 

0.433 

3.14 

81.3 

0.256 

t=0.065" 

56 

0.384 

2.6? 

72.3 

0.271 

0=0. 20 

57 

0.336 

2.16 

63.2 

0.292 

58 

0.288 

1.76 

54.1 

0.324 

Not  Enhanced 

59 

0.240 

1,39 

45.1 

0.368 

60 

0.192 

1.06 

Bsl 

36.1 

0.439 

61 

0.144 

0.75 

27.0 

0.533 

6? 

0.096 

0.47 

75.0 

18.0 

0.779 

63 

0.048 

0.24 

75.0 

9.0 

1.592 

64 

0.125 

0.61 

75.0 

23.5 

0.596 

65 

0.100 

0.47 

75.0 

18.8 

0.718 

66 

0.075 

0.34 

75.0 

14.1 

0.923 

67 

0.050 

0.25 

75.0 

9.4 

1.527 

68 

o.n?5 

0,1? 

75.0 

4.7 

2.933 

18 

TO 

231.3 

0.203 

27  plates 

19 

0.735 

5.48 

75.0 

184.3 

0.204 

0.477 

3.5^ 

75.0 

119.6 

0.310 

d=0.0?0'' 

?1 

0.429 

2.98 

75.0 

107.5 

0.324 

t=0.065" 

?? 

0.381 

2.46 

75.0 

95.5 

0.339 

0=0.20 

?3 

0.333 

2.0 

75.0 

83.5 

0.360 

?4 

0.286 

1.6 

75.0 

71.6 

0.392 

Enhanced 

?5 

0.238 

1.24 

75.0 

59.6 

0.438 

?6 

0.190 

0.92 

75.0 

47.7 

0.508 

?7 

0.143 

0.65 

75.0 

35.7 

0.639 

28 

0.095 

0.41 

75.0 

23.8 

0.907 

?9 

0.048 

0.20 

75.0 

11.9 

1.77 

30 

0.124 

0,52 

75.0 

31.0 

0.678 

31 

0.100 

0.40 

75.0 

24.8 

0.815 

32 

0.074 

0.30 

75.0 

18.6 

1.086 

33 

0.050 

0.20 

75.0 

12.4 

1.630 

34 

0.025 

0.10 

75.0 

6.20 

3.260 

Table  A2.  Continued 
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APC'HpO) 

VluidC^F) 

'^Re 

f 

GAS 

Core 

Description 

mmtyjm 

8.28 

231.2 

056 

Ai  r 

27  Plates 

BiiGS 

5.26 

75.0 

184.3 

0.196 

■am 

3.17 

119.5 

0.279 

d=0.020" 

38 

0.429 

2.66 

75.0 

107.5 

0.290 

t=0.065" 

39 

0.381 

2.21 

75.0 

95.5 

0.305 

a=0.20 

40 

1.79 

75.0 

83.5 

0.323 

41 

0.286 

1.44 

75.0 

71.6 

0.353 

Not  Enhanced  ! 

42 

0.238 

1.12 

75.0 

59.6 

0.3% 

1 

43 

0.190 

0.84 

75.0 

47.7 

0.464 

j 

44 

0.143 

0.59 

75.0 

35.7 

0.580 

i 

45 

0.095 

0.37 

75.0 

23.8 

0.818 

! 

46 

0.048 

0.19 

75.0 

11.9 

1.681 

' 

47 

0.124 

0.48 

75.0 

31.0 

0.626 

i 

48 

0.099 

0.38 

75.0 

24.8 

0.774 

49 

0.074 

0.28 

75.0 

18.6 

1.014 

50 

0.50 

0.19 

75.0 

12.4 

1.548 

51 

0.245 

0.10 

75.0 

6.2 

3.260 

86 

6.32 

tto 

174.3 

rt.l38 

32  Plates 

87 

0.742 

4.14 

77.0 

139.1 

0.141 

88 

0.482 

2.61 

77.0 

90.3 

0.211 

d=0.015" 

89 

0.433 

2.20 

77.0 

81.2 

0.220 

t=0.045" 

90 

0.385 

1.85 

77.0 

72.2 

0.234 

a=0.20 

91 

0.337 

1.53 

77.0 

63.1 

0.253 

92 

0.288 

1.23 

77.0 

54.1 

0.277 

Alternatina 

93 

0.240 

0.97 

77.0 

45.1 

0.315 

Patterns 

94 

0.192 

0.73 

77.0 

36.0 

0.370 

95 

0.144 

0.54 

77.0 

27.0 

0.487 

96 

0.096 

0.34 

77.0 

18.0 

0.690 

1 

97 

0.048 

0.18 

77.0 

9.0 

1.461 

1 

98 

0.125 

0.42 

77.0 

23.5 

0.502 

1 

99 

0.100 

0.34 

77.0 

18.8 

0.635 

100 

0.075 

0.25 

77.0 

14.1 

0.831 

101 

0.050 

0.18 

77.0 

9.4 

1.346 

102 

0.025 

0.10 

77.0 

L  - - 

4.7 

2.991 

Table  A2.  Continued 
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Run 

1 

APC'HpO) 

V]u1d(°F) 

"Re 

' 

f 

GAS 

Core 

Description 

103 

0.926 

TTO 

231.5 

0.174 

Ai  r 

32  Plates 

104 

0.739 

77.0 

184.7 

0.180 

105 

0.480 

77.0 

119.9 

0.256 

d=0.0?0" 

0.432 

1.99 

77.0 

107.9 

0.265 

t=0.045" 

107 

0.383 

1.64 

77.0 

95.8 

0.277 

a=0.20 

108 

0.335 

1.33 

77.0 

83.8 

0.293 

109 

0.287 

1.05 

77.0 

71.8 

0.315 

Al ternating 

110 

0.239 

0.81 

77.0 

59.8 

0.350 

Patterns 

111 

0.191 

0.60 

77.0 

47.9 

0.406 

112 

0.144 

0.41 

77.0 

35.9 

0.493 

113 

0.096 

0.26 

77.0 

23.9 

0.704 

114 

0.048 

0.13 

77.0 

12.0 

1.407 

115 

0.125 

0.34 

77.0 
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g 

APC'H^O) 

Vluid(°F) 

"Re 

f 

GAS 

- r 

Core 

Description 

u] 

70 

170 

0.354 

Air 

29  Plates 

143 
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0.050 

0.13 
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163 
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0.07 
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5.128 

Table  A2,  Concluded 
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APPFNniX  B  -  ACCURACY 


Roth  the  experimental  errors  inherent  to  the  single  blow  transient 
test  method  and  the  inaccuracies  particular  to  the  experimental  set  up 
have  been  considered  [H  to  estimate  the  following  uncertainties: 

Quantity  Maximum  Uncertainty  (percent) 


Np^  Prandtl  Number  ±5.0 

air  Specific  heat  ±  0.5 

.  Dynamic  viscosity  ±  1.0 

C  ®  Specific  heat  ±  4.0 

dP  Hydraulic  diameter  ±  ?.0 

Np  Hole  Reynolds  number  ±  0.0 

T  ^  friction  factor  ±10.0 


Number  of  heat-transfer  units  ±20.0 
Table  Bl.  Maximum  uncertainty  of  physical  constants  and  measurements 

With  reguard  to  equipment  inaccuracies,  some  verification  of  actual  performance 
was  made  where  possible.  Thermocouples  were  checked  and  found  to  give  accuracies  of 
±  0.1  °F.  Three  Fisher-Porter  rotometers  were  used  to  determine  mass  flows.  The 
smallest  meter  was  checked  with  a  wet-test  meter  and  found  to  be  accurate  to  within 
5  percent.  The  estimation  of  maximum  slope  itself  has  an  error  which  varies 
according  to  the  range  of  Reynolds  number  under  consideration.  At  low  Reynolds 
numbers  possible  error  in  reading  maximum  slope  means  a  ±  6  percent 
variation  in  per  plate.  At  higher  Reynolds  numbers  this  reduces  to  ±  3  percent 
By  digitizing  the  time-temperature  profile  through  the  use  of  a  minicomputer  this 
error  could  be  reduced  to  less  than  ±  1  percent.  This  was  done  for  the  eiqh*^ 
electron  beam  drilled  plates. 
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APPENDIX  C  -  MODELINR  THE  PLATE 

An  analytical  model  of  the  plate  was  developed  to  evaluate  the  experimental 
data.  The  Tin  effectiveness  for  a  circular  fin,  the  effective  thermal 
conductivity,  the  friction  factor  vs  hole  Reynolds  number,  and  the  number  of 
heat-transfer  units  per  plate  vs  hole  Reynolds  number  were  determined 
analytically. 

At  first  the  plate  was  modeled  as  a  "matrix,"  similar  to  a  porous  media,  but 
the  correlation  parameters  generated  didn't  describe  the  system  correctly. 

Finally,  the  plate  was  modeled  using  the  staggered  (hexagonal  1  hole  pattern  shown 
in  Figure  Cl.  This  model  adequately  describes  the  heat  exchanger's  behavior. 


t  *  plate  thickness, 
d  *  hole  diameter, 
b  *  base  of  hexagon, 
h  =  height  of  one 

equilateral  triangle 
in  the  hexagon, 
s  =  distance  between 
plates. 


Figure  Cl.  Model  of  staggered  hole  pattern 
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1.  DETERMINE  FIN  EFFECTIVENESS  FOR  A  CIRCULAR  FIN 

The  sinqle-blow  transient-test  method  generates  the  N^^  for  a  plate  with  100% 
fin  effectiveness  (o^)  because  each  plate  is  at  a  constant  temperature  at  any 
given  time  during  the  test.  The  fin  effectiveness  must  be  considered  to 
determine  the  actual  N^^/plate  when  designing  heat  exchangers. 


1 

1 

+ 

\ 

(Cl) 

"tv. 

'^tULp\p 

HP 

The  two  pass  counterflow  heat  exchanger  (see  Figure  2  and  C?)  will  be 
designed  with  a  high  pressure  central  core  and  a  low  pressure  stream  in  the 
annular  section  surrounding  the  center. 


Fiqure  C2.  Two  pass  heat  exchanger 
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Kays  and  London  [31  have  evaluated  and  presented  the  fin  effectiveness  for 
annual  flow  in  Fiqure  2-11.  The  fin  effectiveness  is  also  needed  for  the  circular 
center  section  (hiqh-pressure  stream  in  RCUBE  application!.  The  analysis  follows. 


Figure  C3.  Model  of  circular  fin  (high-pressure  flow  path) 


Start  with  the  differential  equation  describing  temperature  profile  of  the 
fin  shown  in  Figure  C3. 


where 


Let 


d^T 


+  1  dT  +  1 
7?F  TTt 


[q(r)l 


0 


(C21 


T  =  Temperature  at  radial  location  r. 

k  =  Coefficient  of  thermal  conductivity, 

h  =  Convective  heat  transfer  coefficient 

q(r)  =  Absorbed  heat  flux. 

=  hfTg  -  T]. 

T^  =  Temperature  of  the  gas. 

T  =  Tg  +  0  and  transform  eq  (C2)  to 


d^9  .  1  de  .  he 

dr^  r  dr  Ft 


0 


(C3) 


i 

I 
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To  solvp  the  differential  equation,  the  boundary  conditions  must  be 
determined.  Figure  C4  shows  the  temperature  profiles  for  both  Tg  <  T  and  T  >  T^. 


T  >  Tg  Tg  >  T 

Figure  C4.  Fin  temperature  properties 


From  Figure  C4,  the  following  boundary  conditions  can  be  seen. 

At  r  *  0 

—  =  0;  9  is  finite 
dr 

At  r  =  r^ 

0 

First,  tind  ^ 

dr  r 

0 

=  -kf^TT^r  t  ^jr 
ro  0  j  '  0 

dr 

=  -kf^TTr  t)  — |r^  fC4l 

^  dr  " 
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But,  the  following  relation  for  is  based  on  the  definition  of  q  : 

ro  ^ro 


0 

'^ro  ^  ^ 

0 


r 

0  • 

=  -Pirh  S'  erdr 

0  (C5) 

Substituting  equation  (C5l  into  equation  (C4): 


r 

0 

Prrh  s  0rdr 


dr  k(?Trr^t) 


.  H  jO  erdr 

dr  0  kr^t  (Ce) 

Since  the  fin  effectiveness  is  defined  as  the  ratio  of  the  actual  heat 
transferred  to  the  heat  that  would  be  transferred  if  the  entire  fin  was  at  the 
base  ternperature. 


n 

f 


r 

0 

; 

n 


iT-T  ) 
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2TTrdr 


r 

0 

; 

n 


(T  -T  )  ?7rrdr 

0  q 
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Pit/  0  rdr 

0 _ 

irr  2(T  -T  ) 
0  0  g 


2  /  erdr 
0 

r  20 
0  0 


Combine  equations  (C6)  and  (C7): 


2  d6,  kr  t 
-2- 

dr  h 


n  = 


r  20 
0  0 


2k  t  d0| 


ro 


hr  9  dr 
0  0 


Define  as  follows: 


i2  = 


2  = 


hr2 

kt 


hr  2 
0 


kt. 


fC7) 


fC8) 


(C9) 
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Pr  deir 
_ 0  _ '  0 

u  J 

0  0  dr 


To  ^ind  — ,  use  the  chain  rule 
du 

du  dr  du 


and  differentiate  equation  fC9)  to  find  —  : 

du 

2u  du  =  —  2r  dr 
kt 


dr  _  u  kt 
du  r  h 


du  0 


^ir  .y. 

dr  0  r  h 


fClO) 


d9  %  de  (Cll) 

— iir  =  —  —  lU 

dr  0  r^  du  0 

0 
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Substitute  equation  (Cll)  in  equation  (CIO): 
u  . 

_ 0  0  de  u„ 

n  =  ■  ^ 

u  20  r  (iu 

0  0  0 

n,  =  ^|u  fCl?^ 

u  9  du 
0  0 

Now,  find  equation  (C?)  in  terms  of  u  and  0: 

First,  multiply  through  by  r^ 

r^  d^Q  +  rd0  -  r^  h  0  =  0 

dr2  dr  kt  (C13) 


Rearrange  and  differentiate  equation  (C9). 


dr2 
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Use  the  chain  rule  to  Find  —  and  . 

du  dr* 


d*e  ^  ^  ^  ^  /  duV  ^  h_  ^ 

dr*  du  dr*  du*  \dr/  kt  du* 

Now,  equation  (Cl 3)  can  be  rewritten: 

u*  ^  +  u  ^  -  u*0  =  0  fC14l 

du*  du 


Equation  (C14)  is  the  modified  Bessel  equation  (n=0)  and  the  solution 
follows: 

6  =  Cjl^(u)  +  c^K^fu) 

Apply  the  boundary  conditions  to  find  the  constants. 

At  r  =  0,  u  =  0,  9  =  finite 

/ 

Since  =  «>,  c^,  =  0. 

e  =  CjTq(u)  (C15) 


d0  c,  d[T„(u)l 
—  =  ^  ^  =  c  I  (ul 

du  du  ^  ^ 


(C16) 
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Substitute  equations  (C15)  and  (C16)  in  equation  (C12).  Equation  (C17)  is  a 
relationship  for  the  fin  effectiveness  as  a  function  of  the  fin  fhigh-pressure 
center  section)  radius,  the  plate  thickness,  the  heat-transfer  coefficient,  and 
the  thermal  conductivity  of  the  plate. 


n 

f 


where 


2.  EVALUATinN  OF  EFFECTIVE  VALUES  OF  TMFP'  <  CONOUCTIVITY 


^C17) 


The  ef-fective  thermal  conductivity  must  be  used  to  calculate  fin 
effectiveness  and  other  heat-transfer  characteristics.  The  holes  in  the  plate 
cause  a  lower  thermal  conductivity  than  the  corresponding  material  thermal 
conductivity.  The  following  analysis  finds  the  effective  thermal  conductivity  as 
a  function  of  percent  open  area.  Consider  the  flow  of  a  constant  amount  of  heat, 
q,  through  the  distance  2b  defined  by  planes  at  1  and  2  (see  Figure  C?). 


t  =  plate  thickness, 
b  =  one  side  hexagon, 
h  =  height  of  one 

equilateral  triangle 
in  the  hexagon, 
r  =  radius  of  hole, 
q  =  heat  -flow. 


Figure  C5.  Model  of  a  single  hole 
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The  followinq  assumptions  were  made: 

a.  There  is  no  heat  "ffow  across  the  boundaries  defined  by  the  hexaqon  (due  to 
symmetry) . 

b.  The  heat  transfer  by  solid  conduction  is  isotropic  for  the  integrated 
effect  of  the  flow  area. 

The  rate  of  heat  conduction  for  steady  one-dimensional  conduction  follows: 


q  =  k  A(x)  — 
dx 


(C18) 


where 

0  =  Rate  of  heat  conduction, 
k  =  Thermal  conductivity  of  the  material. 
dT/dx  =  Temperature  qradient  at  the  section. 

X  =  Distance  along  the  heat  flow  lines;  defined 
positive  in  the  direction  of  the  heat  flow. 
A(x)  =  Area  the  heat  is  flowing  throuqh  measured 
perpendicular  to  the  heat  flow. 

NOTE:  The  area  is  a  function  of  x. 

Rearrange  equation  (CIS)  and  integrate: 

Ax 

^  ^ 

k  Q  A(x) 

(C19) 

=  1 


AT 

/  dT  = 
0 
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where 

AT  =  The  temperature  chanqe  across  the  section 
from  1  to  2. 

Ax  =  The  distance  along  the  heat  flow  lines  from 
section  1  to  2  ffor  one  cell). 

^  =  Average  value  of  the  reciprocal  of  the  heat 

Afx)|  flow  area. 

To  determine  the  effective  thermal  conductivity, 
let  Ax  =  2h. 

Because  the  holes  are  staggered,  the  area  of  the  two  cells  must  be  considered  when 
Ax  =2b  (see  Figure  C6l. 


o  @ 


Figure  C6.  Model  of  two  holes 
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The  average  area  for  the  coupled  system  of  holes  is  twice  the  area  for  one 
hole  and  the  plate  thickness  is  constant;  therefore,  equation  fC19)  can  he  written 
as  follows: 


^  _  q  1  _  _g _ 1_ 

9b  2kA{x)lgy  2kty(x)|3^ 

(C201 


where 


yfx)  =  The  flow  cross  section  as  a  function  of  x 
for  both  parallel  flow  paths 
(see  Figures  C5  and  C7). 


Equation  (C20)  uses  the  actual  area  and  the  actual  thermal  conductivity.  The 
effective  thermal  conductivity  is  defined  as  the  thermal  conductivity  that  would 
be  needed  to  give  the  same  results  as  equation  fC20)  if  there  were  no  holes.  I'f 
the  effective  thermal  conductivity  and  the  area  assuming  no  holes  are  used  in 
equation  (C19),  the  following  equation  results: 


^  _  _g _ 

2b  3.5tbk' 


(C211 


where 


k'  =  effective  thermal  conductivity 

Next,  find  the  ratio  of  effective  to  actual  thermal  conductivity  by  combining 
equations  (C17)  and  (C21). 

^  _2 _ 

k  3.5b  J, 

y(x)l3y  (C22) 
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Now,  if  [y(x)3av  is  redefined  to  be  the  cross  section  for  a  single  parallel 
tlow  path  instead  of  both  parallel  flow  paths  as  shown  in  Figure  C?,  equation  fCS) 
will  become 


iL  =  4 

k  3.5b  I 

y(x)lav 

(C231 


The  relationship  of  k'/k  as  needed  as  a  function  o^  a  plate's  open  area  ratio 
fol.  To  find  k'/k  as  a  function  of  a,  be  found  as  a  function  of  a. 
First,  consider  a  sixth  of  a  hole  and  the  equilateral  triangle  associated  with  it 
(see  Figure  C7). 


b 


=  One  side  of  the  hexagon  or 
equilateral  triangle. 

=  Radius  of  hole. 

=  Distance  along  the  radius 
from  the  outside  of  the  hole 
to  the  edge  of  the  hexagon. 

=  Angle  to  the  radius  measured 
from  the  height  of  the 
equilateral  triangle. 


Figure  C7. 


Model  of  a  sixth  of  a  hole 


The  hole  geometry  determines  cos 6. 


COS0 


b 


r  +  yix) 


rc?4l 
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y(x) 


-JP ,  b 


COS0 


r 


JL_  1  (C25) 

y(x) 

COS0 


Now,  r/b  must  be  found  as  a  function  of  a.  a  1s  the  ratio  of  plate  open  area 
to  plate  frontal  area  and  Is  also  the  ratio  of  the  area  of  one  hole  to  the  hexagon 
surrounding  that  hole. 


6[fl/2)b3[(^/2)b] 


2irr2 


b2 


r/b  =  0.90<)4 


(C26) 


Substituting  equation  (C26)  In  equation  (C25)  results  In  the  following 
relationship: 


=  (&.  .  0.9094/o^)"^ 
y(x)  COS0 
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b 

For  anv  given  a,  can  be  found  by  taking  the  average  value  for 

O^0<3O. 


b 


3'^>'^lav 


6 

_  j  (J?H  .  0.9094^)"^  d0 

Tf  0  COS0 


(C28) 


Solving  equation  (C?8)  for  various  values  of  o  and  plugging  the  result  into 
equation  (C22l  will  generate  k'/k  (a).  Figure  C8  shows  k'/k  plotted  vs  o. 


Figure  C8.  Plot  of  effective  thermal  conductivity  ratio  vs  plate  ofen 
area  ratio 

As  the  open  area  approaches  zero,  k'/k  approaches  unity  as  it  should;  it 
doesn't  go  exactly  to  one  because  of  the  approximations  involved  in  the  system 
analysis. 
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This  analysis  was  confirmed  experimentally  by  using  the  following  analogue 
between  thermal  conductivity  and  electricity: 

k'  ^  ^solid  (C29) 

k 

where: 

V  =  Voltaoe  across  solid 
solid 

material  for  a  given  current. 

V  =  Voltage  across  perforated 

material  for  the  same  current. 

3.  DETERMINE  FRICTION  FACTOR  vs  HOLE  REYNOLDS  NUMBER 

A  single  plate  will  be  considered  to  find  the  friction  factor  for 
the  entire  core.  Figure  C9  shows  the  plate  to  be  considered  and  the 
corresponding  pressure  drops. 


Figure  C9.  Plot  of  pressure  drop  across  one  plate 
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The  change  in  pressure  across  the  plate  is  a  combination  ot  entrance  static 
pressure  losses,  friction  and  acceleration  losses,  and  exit  static  pressure 
recovery  for  each  hole.  Equation  (C30)  quantifies  these  pressure  changes  as 
presented  by  Kays  &  London  [3]. 


AP  Vj,  A  v^ 

_=  -  [fl+K-a2)  +  2{ - 1)  +  - -  (l-K-a^— ] 

c  ^ 


1 


1 


'1 


'1 


entrance 

loss 


acceleration 

loss 


friction 

loss 


exit  (C30) 
reqa i n 


where 

AP  =  Pressure  drop  across  the  plate. 

Pj  =  Pressure. 

G  =  Flow  stream  mass  velocity. 

V  =  Specific  volume. 

q^  =  Proportionality  factor  in  Newton’s  second  law. 

=  Entrance  coefficient  as  defined  by  Kays  ft 
London  [3]. 

Kg  =  Exit  coefficient  as  defined  by  Kays  A  London  TS 

a  =  Ratio  of  open  area  to  frontal  area. 

A  =  Total  heat-transfer  area  of  one  hole. 

A  =  Minimum  free  flow  area  of  one  hole, 
c 

f  =  Mean  friction  factor  in  holes. 

Subscripts 

1  =  Value  at  plane  1  (Figure  C9l. 

2  =  Value  at  plane  ?  (Figure  C9). 
m  =  Average  value. 
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The  friction  factor  (f)  is  defined  by  the  following  relationship.  This 
observed  friction  factor  includes  entrance  and  exit  losses. 


AP  =  f‘  ^ 

L 

’"h 

where 

P 

=»  fluid  density. 

L 

=  Length  of  test  section. 

"h 

=  A^L/A  =  hydraulic  radiu 

Since  only  one  plate  is  being  considered,  the  total  length  of  the  exchanger 
is  .iust  the  thickness  of  a  plate  (L  =  t). 


A'P 


G^v  A 

m  f ' 


AP  = 


f  A  ''m 


(C311 


?g. 


A  ''1 


Equating  equations  (C30)  and  (C31)  results  in  a  relationship  for  the  observed 
friction  factor  per  plate: 

-  1)  -  ^(A/A^)^„,^(v^V;)  - 


(A/A  )  ,  .  (v  /v, ) 
'  c 'plate  '  m'  1' 


(C3?) 


The  hydraulic  diameter  is  related  to  the  hydraulic  radius  and  A^/A: 

n.  4r.  4  A 

h  _ _ n  _  _c 

L  L  A 
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Rearranging  this  results  in  the  following  relationship  for  A/A^: 

A  = 

A  n.  fC33) 

c  h 

For  any  circular  cylinder,  the  diameter  of  the  cylinder  is  the  hydraulic 
diameter,  and  L  is  the  plate  thickness  when  only  one  plate  is  considered. 
Substituting  these  values  in  equation  (C33)  and  substituting  equation  fC33)  in 
equation  fC32)  results  in  the  following  equation  for  the  observed  friction  factor 
for  one  plate: 

[(l+K^-a2)  +  2(v^/vj  -  11  -  4f(t/d)(v^/vj)  -  (1-0^-  K^){vp/Vjl] 

f  = 


(Vl' 


fC34) 


To  determine  the  friction  factor  as  a  function  o'*"  hole  Reynolds  number,  the 
various  components  of  equation  fC34)  must  be  determined.  First,  will 
be  determined  using  the  plate  model. 


where 


A/A 

c 


A/A 


fr 


A^^  =  Total  frontal  area. 

A  =  Total  heat-transfer  surface  area. 
D  =  Diameter  of  test  section. 


Let  =  Area  of  hexagon  around  the  hole. 

Ao  =  Area  of  the  hole. 

Total  surface  area  of  connected  voids 

a  =  — - 

Total  volume  of  matrix 
=  area  density  f see  Figure  Cl). 
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The  area  density  definition  is  used  to  find  the  area  density  as  a  function  of 
hole  geometry  and  plate  open  area. 


■  X  0-' 


a  = 


2 


A^  4Trrt 
0 


t  A 


b^t 


=  fl-a)  + 


Let  A'  and  A'^  be  the  areas  corresponding  to  a  single  hole. 


A'  =  aA^t  =  at  o/q 


t  r  o 


A  '  =  TTr^ 
c 


A  ,  A^  ,  opfl-o)  ^  t 


fC35) 


Next,  ^ind  v^/v^ 


For  a  pt  rfect  gas  v=  ^ 

P 


''I  ^  ^2 
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Let 


aP  =  The  pressure  drop  across  a  single  plate, 
n  =  Number  of  plates. 


Then 


EAP 

n 


=  Average  pressure  drop/plate. 


For  a  given  plate  the  pressures  can  be  related  to  the  average  pressure  drop 
per  plate. 


P 


O 


eap 

n 


P,/P,  =  1  - 


?:ap 

np. 


Now,  for  any  reasonable  design,  the  pressure  drop  across  the  whole  heat 
exchanger  will  be  less  than  five  percent,  and  there  will  be  more  than  ten  plates 
so,  for  a  reasonable  design  Pp/P|  =  1  because 


EAP 


nP 


1 


<  0.005 


Since 


Pp/P 


1 


s  1 
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The  same  reasoninq  is  used  to  find 


AT  =  Temperature  drop  across  an  exchanger  having  n  plates. 


=  1  ±  ^ 


If  the  stream  is  being  heated  T^/T^  s  1  and  if  it  is  being  cooled  Tj,/Tj  5  1 


so,  Vj,/v^  can  have  two  valves: 


Vp/Vl  =  1  ±  ^ 


(C36) 


v^  is  the  mean  specific  volume.  v„/v,  can  be  found  the  same  way  v^/v, 
m  ml  ?  1 


was  found. 


V  /V.  =  ^T^/T.1  fP./P^l 

ml  ml  1  m 


Pl/Pn,  =  ■' 
1  m 


T  =  T,  ± 
m  1 


V"i  .  1  ± 


V'-I  '  '  *- 


(C371 
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Recause  of  temperature  variations  alonq  the  length  of  the  exr.hanqer,  the 
value  of  v^/Vj,  and  therefore  f,  will  vary  -from  plate  to  plate.  Average  values 
of  Vj>/Vj  and  v^/v^  will  be  used.  Vfhen  these  values  are  applied  to  one  plate  and 
multiplied  by  the  number  of  plates  in  the  exchanger,  the  resulting  total  pressure 
drop  will  be  the  pressure  drop  measured  over  the  whole  exchanger. 


AT 

1 

^0? 

AT 

—  = 

— 

/ 

—  dT 

nTj |av 

AT 

^ol 

AT 

1 

"o? 

=  — 

In 

nT^ 

|av  n 

^01 

(C3R) 


where 


TOp  =  Temperature  exiting  hea+  exchanger. 


To 


1 

AT 


Temperature  entering  heat  exchanger. 

Jop  -  To^ 

Temperature  at  the  inlet  of  any  given  plate 


The  average  values  of  v^/v^  and  v^/v^  follow: 


lav  =  1  in 


To, 


To, 


fC39) 


V  /v,  |av  =  1  +  -  In 
m  1  '  - 


?n 


To 


1 


(C40) 


Next,  the  relationship  between  the  calculated  friction  factor  and  the 
theoretical  friction  factor  must  be  found. 
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Let 

t  =  Measured  friction  factor, 
e 

f  =  Theoretical  observed  friction  factor. 


The  measured  friction  factor  calculated  from  the  data  is  defined  as  follows: 


•f 

e 


n.  AP 

h  m  c 
2nt 


^C41) 


Also,  the  followinq  theoretical  relationship  exists  for  AP  over  one  plate: 


AP 


n2v 


m 


f 


A 

A 


c 


fC4?) 


Substitute  equation  (C4P)  in  equation  (C41).  Equation  (C42)  is  the  pressure 
change  over  one  plate  so  n  =  1  for  equation  (C411.  The  Reynolds  number  is  based 
on  the  hole  diameter  so  =  d. 


2t  2g^  A^  v^G2 

fp  ^  _d  f '  A 
4t  A 

c 


(C43) 


Substituting  equation  fC34)  into  equation  fC43l  results  in  the  followinq 
expression  to  theoretically  estimate  the  measured  friction  factor: 


=  rd/4t)rVj/v^)r(l+K^-a2l  +  ?^Vj,/Vj  -1)  + 

4f(t/d)(Vpj/Vj)  -  (l-a^-KpKvp/Vj)!  (C44) 
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The  test,  to  detemine  the  friction  factor  was  performed  at  steady  state  room 
temperature.  To  compare  the  theoretical  values  with  the  experimental  values,  the 
following  relationships  will  be  used: 


For  these  conditions,  equation  (C441  becomes 


*  K^)  +  f  . 

4t 


The  friction  factor  per  plate  is  a  function  of  plate  geometry  and  entrance, 
exit,  and  flow  conditions  in  the  hole.  Kays  and  London  present  and  in 

Figure  5-P  T.?]  and  f  in  Fioure  6-2?  [31.  In  both  Tables  the  Reynolds  number 
is  based  on  the  hole  diameter  so,  the  hole  Reynolds  number  (Np^)  is  the  same  as 


Table  Cl  presents  the  friction  factors  per  plate  for  both  punched  and  etched 
plates  calculated  using  equation  IC44l.  Table  C2  presents  the  theoretical 
friction  factors  for  the  ERn  plates. 
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Plate  Description 

Nd 

Re 

NR^(d/t) 

•k 

4(t/rtVn,^ 

^c 

H 

fl 

Punched  Plates 

10 

8.571 

16 

0.467 

0.97 

0.45 

1.9 

d  =  0.018" 

20 

17.14 

18 

0.233 

0.97 

0.46 

1.2 

t  =  n.n?i" 

50 

42.86 

21 

0.093 

0.71 

0  =  n.?36 

Kg 

85.71 

25 

0.047 

0.78 

0.47 

0.52 

■ 

171.43 

32 

0.023 

0.75 

0.48 

0.42 

Ftched  Plates 

■D 

15.83 

17 

0.253 

1.2 

0.76 

2.5 

d  =  0.0285" 

Hi 

31.67 

20 

0.126 

iga 

0.76 

1.7 

t  =  0.018" 

50 

79.17 

24 

0.051 

EM 

0.78 

1.1 

0  =  0.09 

100 

158.3 

30 

0.025 

EM 

0.79 

0.93 

200 

316.7 

38 

0.013 

0.73 

0.79 

0.79 

Ftched  Plates 

10 

15.95 

18 

0.251 

0.98 

0.43 

2.4 

d  =  0.0295" 

90 

31. 8Q 

21 

0.125 

0.97 

0.43 

1.6 

t  =  0.0185" 

79.73 

25 

0.050 

0.80 

0.44 

0.99 

0  =  0.245 

159.46 

31 

0.025 

0.70 

0.48 

0.78 

H 

318.92 

38 

0.013 

0.65 

0.48 

0.64 

*  Read  from  Kays 

and  London 

,  Table  6 

-22  rai. 

**  Read  from  Kays 

and  London 

,  Table  5 

-2  r3i 

• 

Table  Cl.  Theoretical  friction  factor  for  punched  and  etched 
plates 
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Table  C2.  Theoretical  friction  factor  for  EBD  plates  20%  open 
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Table  C^.  Concluded 

4.  nFTFRMTNF  N^^^/Plate  VS  HOLE  REYNOLDS  NUMBER 

Once  the  N^^/plate  is  found  as  a  function  of  Reynolds  number,  the  Colburn 
modules  can  be  calculated  usinq  equation  (1).  The  N^^/plate  was  used  to  compare 
the  experimental  and  theoretical  results. 

N  ^  HA  fC45l 

tu  iTlc„ 

P 

wherp 

II  =  Overall  thermal  conductance. 

A  =  Total  heat-transfer  area. 

rfi  =  Mass  flow  rate. 

c  =  Specific  heat  of  fluid. 

P 
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To  find  HA,  divide  the  total  heat-transfer  area  into  three  parts:  the  front  ^ace, 
the  cylindrical  hole,  and  the  back  face.  Once  the  area  averaged  heat-transfer 
coefficient  is  found  for  each  area,  UA  can  be  found  from  the  following 
relationship 

DA  =  hjA^  +  hpA^  +  h^A^ 

Where  h  =  Area  averaged  heat-transfer  coefficient. 

Subscripts 

1  =  Front  face. 

7  -  Cylindrical  hole. 

3  =  Rack  face. 

UA  can  be  approximated  by  assuming  the  front  and  back  surfaces  have  the 
same  heat-transfer  characteristics. 


UA  =  2hjAj  +  hjjAj, 


fCAF) 


A^  is  the  solid  surface  of  the  plate,  which  is  the  frontal  area  times  the  percent 
sol  id  area . 


Aj  =  n  -  a)  fC47) 

The  total  heat-transfer  area.  A,  can  be  found  by  manipulating  equation 
(C35). 

A  =  +  ii 

A^  a  r  (C35] 

A  =  *  ^1A 

a  d  C  A 

c 


A  =  2A^^ri_Q  +  gtol 

d  (C48) 
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*56  found  from  the  following  relationship: 
A  =  ?A^  +  Aj, 


(r,49) 


^2"  ^ 

-  2A^^ri“0-:^1  -  ?A^^(1-ct) 

d 


4A^^ta 


d 


(C5nl 


Next,  tind  an  expression  for  itic 


P 


as  a  function  of  hole  Reynolds  number. 


'"^p  pVn  iti  A^ 
•  -  ^Re^^^fr^p 


where 


N 


Re 


u 

p 


V 

n 

A 

I 

a 


Hole  Reynolds  number. 

Viscosity. 

Fluid  density. 

Fluid  velocity  upstream  of  ■f'est  section. 
Diameter  of  tubing  upstream  of  test  section. 
Open  area  of  plate. 

Ratio  of  open  area  to  plate  frontal  area. 
Frontal  area  of  plate. 

Hole  diameter. 


(C51) 
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Mow,  substitute  equations  fC47),  (C41),  and  (C51)  in  equation  fC46). 


(hjlPA^^d-a)]  +  h2(A^^at/d)}d 

N_  uoA^  c 
Re^  tr  p 


^  1  r2hjd(^)  +  4hpt] 

^tu  ^  (C52) 

Now  it's  necessary  to  find  and  h^.  b^  can  be  found  as  a  function  of  the 

upstream  mass  velocity  ^before  the  flow  enters  the  plate),  R',  and  the  equivalent 
cylinder  diameter,  d'.  Due  to  continuity,  the  upstream  mass  velocity  is  related 
to  the  mass  velocitv  in  the  core  as  follows: 

R'  =  Ga  fC53l 


The  eouivalent  cylinder  diameter  can  be  approximated  by  the  following  eouation; 


d'  =  2yfx) 


av 


fC54) 


Note  that  y(x  is  shown  on  Figure  C6. 


d’  2y(x) 


av 


yd) 


av 


fC55) 


Substitute  equation  (C24)  fnry^x)|  in  equation  fC55l,  then  replace 

B  V 

r/b  with  equation  rC26); 


d'  _  rJ^/  ? 

1  i 

d  0.Q094o  cose 
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Integrate  over  one  12 


th 


of  a  hole  to  find  d'/d: 


d' 

d 


6 


IT 


it/6 


/  ( 
0 


n.95? 


1 

a  cose 


-  Dde 


fC56l 


A  relationship  for  d'/d  as  a  function  of  a  can  be  approximated  by  correlating 
values  evaluated  using  equation  (C56).  Equation  (C57)  is  accurate  to  10%  of  d'/d 
obtained  by  numerically  integrating  equation  {C56)  for  a  S  0.3. 


d'/d  =  .358 


{C57) 


A  relationship  for  h^  (C58)  is  presented  in  McAdams'  "Heat  Transmission. "[4'' 
Although  this  relationship  is  based  on  data  from  air,  it  can  be  used  for  the 
helium  heat  exchangers  because  the  Prandtl  numbers  are  approximately  equal.  The 
relationship  is  also  an  approximation  because  equation  (C58l  is  for  flow  normal  to 
a  single  cylinder.  This  was  used  to  approximate  the  f^ow  normal  to  a  flat  plate 
because  no  reference  could  be  found  for  low  Reynolds  number  flow  normal  to  a  flat 
plate.  Or  Fowle  Til  believes  this  is  a  conservative  estimate  —  it  will  predict 
lower  than  actual  N^^/plate. 


hjd' 

k 


1  i P I  0 . 52 
0.32  +  0.43(2-2-1 

V 


for  0.1  <  s  1000 


u 


(C58) 


Note  the  following  relationship  for  the  Reynolds  number: 

^  =  ^Re  fC59l 

V 
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Substitiitina  equations  (C53),  (C57),  and  (C59)  in  equation  fC58)  results  in  the 
followinq  expression  for  h^: 

hj  =  (k/d)  2.793  cr^-'^68  ^  0. 25?Nj^p'^* (060) 

Two  equations  (C61  and  C62)  to  calculate  h.,  are  shown  in  "Heat,  Mass  and 
Momentum  Transfer"  by  Robsenow  and  Choi  TS''.  These  functions  are  for 
heat  transfer  in  a  circular  cylinder  and  are  based  on  the  following  assumptions: 

1  s  Np  <  inno 

"^PR  =  0.7 

Uniform  wall  temperature 

The  first  expression  for  h^  assumes  a  developing  velocity  profile  in  laminar  flow. 

hp  =  k/d  r3.66  +  0.073  d/t  Np^  ^ 

1  +  0.012  (d/t  ^Re)°‘® 

The  second  expression  for  hp  assumes  a  parabolic  velocity  profile  in  laminar  flow. 

h^  -  k/d  (3.66  +  0.047  (d/t  Np^) 

1  +  0.032  (d/t 

Substituting  equations  (C60),  fC61),  and  (C62l  in  equation  (C.5?)  and  recognizing 
+  hat 


uc. 


^Pr  =  — H  =  0.7 

k 


results  in  equations  (C63)  and  (C67). 
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FOR  A  DEVELOPING  VELOCITY  PROFILE  IN  LAMINAR  FLOW 

N^^/Plate  =  fC63) 

-0.232  0.52  O.Pl 

(1-ct)c  12.554+2.011  N^^  a  )  20.9  0.417 

-  + - + - 

■^Re  ^Re  l+.012f'^/t^Re)^-® 

For  the  punched  plate  -  0.18",  t  =  0.021", a  =  0.236) 

N^jj/plate  =  27.1  +  1.803  (C64) 

1  +  0.011  N^g  0.8 

For  the  etched  plates  (d  =  0.0285",  t  =  0.0185",  a  =  0.09) 

N^  /plate  =  17.63  N.  +  2.391  ’^Re'^*^^  +  fC65) 

tu  Re  _ 

1  +  O.Ol/Npg*^-^ 

For  the  etched  platesid  =  0.0295",  t  =  0.018",  a  =  0.245) 

N^jj/plate  =  15.42  +  1.775  Np^  rC66) 

1  +  O.OlBMpp^'^ 
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FOR  A  PARABOLIC  VELOCITY  PROFILE  IN  LAMINAR  FLOW 


tu 


/plate  = 


(C67) 


-0.232 


0.52  0.121 


{l-a)o  (2.554+2.011  a  )  20.9  0.269 


'Re 


^Re  ^^/t^Re^ 


0.667 


For  the  punched  plate 

N^^/plate  =  27.1  +  1.803  ^  _  (C68^ 

1  +  0.029  No 
Re 

For  the  etched  plate  (9^  open  area) 

N^^/plate  =  17.63  +  2.391  ^  _  (C69) 

1  +  0.043No 
Re 


For  the  etched  plate  (24.5%  open  area) 

N^^/plat.e  =  15.4?  N^^'^  +  1.775  _  (C70) 

1  +  0.044N_ 

Re 

Table  C3  presents  the  predicted  values  of  N^^/plate  vs  hole  Reynolds  number  for 
the  developing  and  parabolic  velocity  profiles  in  the  punched  and  etched  plates. 
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Nt„/ptate 

Developing  Velocity  Profile 

Parabol ic 

Velocity 

Profile 

Punched 

Etched 

Punched 

Etched 

a  =0.236 

a  =0.09 

o'  =0.245 

d  =0.236 

q  =0.09 

Nn 

Re 

. 

m 

in 

3.697 

2.931 

2.505 

1 

♦ 

3.544 

2.779 

?n 

2.155 

1.801 

1.541 

2.005 

1.654 

50 

1.151 

1.019 

0.876 

1.011 

0.818 

0.748 

0.759 

0.687 

0.593 

0.634 

0.578 

0.486 

0.514 

0.468 

0.403 

0.412 

0.385 

0.323 

0.307 

0.277 

0.238 

0.241 

n.2?9 

0.191 

Table  C3.  N^^/plate  vs  hole  Reynolds  number  for  the  punched 
and  etched  plates 

The  parabolic  velocity  profile  generated  N^^^/plate  values  closest  to  the 
experimental  results.  Consequently,  a  parabolic  velocity  profile  was  used  to 
compare  the  theoretical  and  experimental  results.  Table  C3  presents  the  N^^^/plate 
vs  Npp  calculated  using  equations  (C56),  (C52),  (C58),  and  (C62).  Instead  of 
using  eo'iation  (C57)  whic^’  estimates  d'/d  as  a  function  o^  a,  equation  (C56)  was 
numerically  integrated  to  find  values  for  d'/d.  The  theoretical  results  were 
correlated  and  b  and  m  were  found  to  make  the  following  equation  fit  the 
theoretical  results: 


N^^/pTate  = 
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Type  of  Plate 

d/t 

0 

\e 

Nj^/plate 

b 

m 

Electron 

0.154 

o.?o 

10 

14.785 

116.63 

-0.907 

Beam 

20 

7.654 

Drilled 

30 

5.253 

40 

4.043 

50 

3.311 

60 

2.820 

70 

2.467 

80 

2.201 

90 

1.992 

100 

1.824 

0.?2? 

0.20 

10 

10.600 

77.78 

-0.878 

20 

5.558 

30 

3.853 

40 

2.991 

50 

2.467 

60 

2.115 

70 

1.861 

80 

1.669 

90 

1.518 

100 

1.396 

0.231 

0.20 

10 

10.251 

74.62 

-0.875 

20 

5.383 

30 

3.736 

40 

2.903 

50 

2.397 

60 

2.056 

70 

1.810 

80 

1.624 

90 

1.478 

100 

1.360 

Table  C4.  Theoretical  N^^/plate  assuming  a  parabolic  velocity  profile 
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Type  of  Plate 

d/t 

0 

^Re 

N^u/P^3te 

b 

m 

n.308 

0.20 

10 

7,982 

54.39 

-0.848 

20 

4.246 

30 

2.976 

40 

2.330 

1 

50 

1.937 

fiO 

1.671 

70 

1.479 

80 

1.333 

90 

1.218 

100 

1.125 

Electron 

0.333 

0.20 

10 

7.458 

49.86 

-0.840  1 

Ream 

20 

3.983 

Dri 1 1 ed 

30 

2.800 

40 

2.198 

50 

1.830 

fiO 

1.582 

70 

1.407 

80 

1.265 

1 

1 

90 

1.157 

100 

1.070 

1 

0.^44 

0.20 

10 

5.886 

36.62 

-0.810 

70 

3.194 

30 

2.271 

1 

40 

1.708 

50 

1.509 

fiO 

1.312 

1 

70 

1.169 

1 

! 

80 

1.060 

90 

0.973 

1 _ 

100 

0.903 

Table  C4.  Continued 
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Punched 


d/t  a 

Re 

Nt„/plate 

b 

n 

0.50n  0.20 

10 

5.362 

32.38 

-0.798 

20 

2.930 

30 

2.094 

40 

1 .665 

50 

1.401 

60 

1.221 

70 

1.090 

80 

0.990 

90 

0.911 

100 

0.846 

0.750  n.po 

10 

3.961 

21.56 

-0.754 

20 

2.224 

30 

1.618 

40 

1.304 

50 

1.109 

60 

0.975 

70 

0.876 

80 

0.801 

90 

0.740 

100 

0.691 

0.857  0.236 

10 

3.567 

18.87 

-0.742 

20 

2.016 

30 

1.474 

40 

1.192 

50 

1.017 

60 

0.896 

70 

0.807 

80 

0.738 

90 

0.683 

100 

0.638 

Table  4.  Continued 
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d/t 

a 

Nn 

Re 

Nj/plate 

b 

Etched 

1 . 583 

0.09 

10 

2.725 

12.41 

-0.674 

20 

1.624 

30 

1.2?4 

40 

1.010 

50 

0.874 

60 

0.779 

70 

0.708 

80 

0.652 

90 

0.607 

100 

0.570 

1.595 

0.245 

10 

2.410 

11.15 

-0.682 

20 

1.426 

30 

1.071 

AO 

0.882 

50 

0.762 

60 

0.679 

70 

0.616 

80 

0.567 

QO 

0.528 

100 

0.495 

Table  C4.  Concluded 
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APPFNniX  n  -  INSULATION  REQUTRFMENTS 

When  a  heat  leak  was  discovered  in  the  test  section,  the  insulation  had  to 
be  chanqed  to  reduce  the  heat  leak.  Althouqh  a  steady-state  analysis  had  shown 
the  insulation  was  adequate,  a,  transient  analysis  showed  that  more  insulation  was 
needed.  A  comparison  ot  the  steady  state  and  transient  analyses  tpllow. 

1 .  STFAOY  STATE 


1.  =  Lenqth  of  test  section. 

r.  =  Radius  of  test  section. 

r  =  Radius  of  test  section  and  insulation. 

0 

h  =  Natural  convective  heat-transfer  coefficient. 

n 

k  =  Thermal  conductivity  of  insulation. 

T^.  =  Temperature  inside  cylinder. 

T  =  Temperature  outside  cylinder. 

a 

Fioure  01.  Schematic  of  test  section 


The  heat  leak  through  the  insulation  can  be  calculated  from  the  following 
equation: 


q  =  2TTr.L(T.  -  Tg)/U 


(nil 
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where 


U  =  (l/hjj)(r^./r^)  +  (l/k)r^.  1n(r^/r^.) 

The  foUowinq  properties  were  used  to  generate  a  curve  of  g  vs  insulation 
thickness , 


r.  =  l.f)  in 
L  =  4.?  in 
h^  =  1  BTII/hr  ft2  F 
AT  =  /T.  -  T  )  = 

k  =  0.n?5  BTU/ft  hr  F  ^polyurethane  foam) 


Applying  these  conditions  to  equation  (Dl)  results  in  the  relationship  shown  in 
Fiaure  D? . 


Insulation  Thickness  (r  -  r.)  in 

o  1 

Figure  0'^.  Pint  of  heat  Teak  rate  vs  insulation  thickness 
(steady-state  analysis) 
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If  h  r./k  >  1  adding  insulafinn  always  d^crpasps  thp  hpat,  Irak.  I-f  h,r./k 
a  1  a  1 

\  therp  is  an  arpa  whpre  adding  insula  +  ion  incrpases  fhp  hpat.  Ipak.  In  this  casp, 
thp  stpadv-state  analysis  shows  that  adding  insulation  would  reduop  thp  hpat  loak 
and  that  1.5"  of  insulation  was  sufficient. 

2.  TRANSTFNT  THERMAL  BEHAVIOR  OF  TMSHLATinN 


k  =  Thermal  conductivity  of 
insulation . 

P  =  Insulation  density, 

c  -  Specific  heat  o^  insulation. 

a  =  k/pc  =  thermal  defussivity, 

s  =  Radius  of  test  section, 

I,  =  Length  of  test  section. 

n  =  Total  hpa+  loss  in  time  t. 
s 

0^=  Temperature  difference 
between  heated  core  and 
ambient  temppra^’urp. 

T=  Time,  measured  from  test 
start. 


Figure  03.  Model  of  the  insulated  test  section 

The  test  section  is  modeled  by  a  circular  cy''indpr  surrounded  bv  ar, 
infinitely  thick  laver  of  insulation  ^Figure  03).  The  temperature  at  the  surface 
of  the  pipe  is  instantaneously  increased  when  thp  test  starts  (t=01.  tL' 

dimensionless  group  aO^/kLs^O^  is  tabulated  as  a  function  of  at/s^  by  .lakob  Ifi}: 
conspguentl  V,  for  a  given  time  (t.)  and  insulator  qualities  (a),  the  total  hea - 
>  loss  can  be  found. 

•  The  test  section  had  the  following  properties: 


s  =  n.l25  ft 
®s  " 

L  =  4.2/1?  ft 
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Table  DI  presents  the  three  insulation  cases  that  were  evaluated. 


_ _ 1 

Case 

1 

Case 

3 

Case 

3 

k 

fRTII/hr  ft  F] 

0.03 

0.02 

0.025 

p 

flbm/-Ft3i 

6 

2 

2.34 

c 

fpTII/lbm  Fl 

0.33 

0.33 

0.33 

T  fserl 

0 

0 

0 

ai 

s 

ax 

s 

CtT 

s 

S2 

(Bill) 

S2 

fBTK) 

S? 

(BTII) 

?0 

0.00533 

0.116 

0.011 

0.055 

0.012 

0.0667 

AO 

0.011 

0.166 

0.021 

0.079 

0.023 

0.0952 

60 

0.016 

0.204 

0.03'’ 

0.097 

0.035 

0.1177 

80 

0.021 

0.238 

0.043 

0.114 

O.046 

0.1383 

100 

0.027 

0.^64 

0.053 

0.130 

0.058 

0.1567 

200 

0.053 

0.389 

0.107 

0.190 

0.115 

0.2290 

300 

0.080 

0.483 

0.160 

0.237 

0.173 

0.2883 

400 

O.107 

0.571 

0.213 

0.28? 

0.230 

0.3396 

800 

0.133 

0.642 

0.267 

0.318 

0.288 

0.3855 

■^ablp  ni.  Comparison  of  di'^terent  insulation  reouirements 
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Fiqurp  D^l  compares  insulation  with  various  kpc  produces.  The  analysis 
indicates  that  total  heat  loss  decreases  as  the  kpc  product  decreases.  Also,  no 
insulation  causes  a  smaller  heat  leak  than  an  infinite  amount  of  insulation  for 
short  runs.  Durinq  most  of  the  test  runs  the  core  reached  steady  state  in  90  to 
660  seconds;  v!nly  one  run  was  under  90  seconds.  The  modified  apparatus  used  3"  of 
insulation  with  a  kpc  product  of  0.019  BTU^/hr  ft^  F  fcase  31. 
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Fiqure  D4.  Transient  analysis  of  heat  loss  for  an  infinite  thickness 
insulation 
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